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Scilab  numbering  policy  used  in  this  document  and  the  relation  to  the 
above  book. 

Exa  Example  (Solved  example) 

Eqn  Equation  (Particular  equation  of  the  above  book) 

AP  Appendix  to  Example  (Scilab  Code  that  is  an  Appednix  to  a  particular 
Example  of  the  above  book) 

For  example,  Exa  3.51  means  solved  example  3.51  of  this  book.  Sec  2.3  means 
a  scilab  code  whose  theory  is  explained  in  Section  2.3  of  the  book. 
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Chapter  1 

Elementary  Concepts  and 
Definitions 


Scilab  code  Exa  l.A.l  Additional  Solved  Problem  1 


t=0 . 1 
R  =  800 

i=50*exp ( - 10*  t ) / 1000 
v  =  50*  cos  (25*t ) 
v_i  =  10*t  ~ (2 . 5) 


Pl=i*i*R 

P2=v*v/R 

P3=v_i 

disp (P3 ,P2 ,P1) 


Scilab  code  Exa  1.A.2  Additional  Solved  Problem  2 
v2  =5 

i2=v2 / 1000 
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vl=-i2/5 
vs  =  vl 

disp (vs ) 


Scilab  code  Exa  1.A.3  Additional  Solved  Problem  3 
L=20/ 1000 

///////////////////////////// 

disp  (”  Parti”  ) 

function  i  =  i (t) 

i  =  8*  exp  (  - 100*  t )  ; 
endf unction 

t=10/1000 

v=  L* ( derivative (i  ,  t)  ) 
d  i  s  p  ( v  ,  ” V=” ) 

///////////////////////////// Part2 
disp  (”  Part2 ”  ) 
t=0 . 1 


function  v=v(t) 

v  =  6*  exp ( -12*t ) 
endf unction 


iO  =8 


i2=l/L*intg (0 , t , v) +i0 
di sp ( i2 ) 


disp  (”  Part3 ’ 
function  i3 


) 

=  i3(t) 
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i3  =  10  - 10*  exp  ( -50*  t )  ; 
endf unction 

t=40/1000 

v3=L*(derivative ( i 3  ,t)  ) 

P40  =  v3  * i3 (t) 
di sp ( P40 ) 

t=50/1000 

P50  =  l/2*L*(i3  (50)  ~2) 
disp(P50) 


Scilab  code  Exa  1.A.4  Additional  Solved  Problem  4 

function  i  =  i (t) 

i  =  9*sin(°/0pi*t)  / 1000  ; 
endf unction 

t=200/1000 
CV= intg (0 , t , i ) 

Wc=300E-6 
C=(CV~2) / (2*Wc) 
disp (C) 


Scilab  code  Exa  1.A.5  Additional  Solved  Problem  5 


L=25E-3 


T=5/ 1000 

I=l/L*integrate(  ’25*008(500*1)  ’  ,  ’ t  ’  ,0,T) 
P=25*cos (500*T)*I 

1 1  =  °/„pi /2/500  /////power  absorbed  zero 
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t2  =  7„pi/500 


// III  energy  stored  zero 


di sp ( I , P , t 1 , t2 ) 


Scilab  code  Exa  1.A.6  Additional  Solved  Problem  6 

eO  =8 . 85E-12 
er  =  1 

A  =  °/0pi*(l/200)  “ 2 
d=150/ 1000000 
C=er*e0*A/ d 

E=l/ 1000 
v=  sqrt (2*E/C) 
disp (v) 

E=2/ 1000000 
v  =  100 

er=2*E/ (v~2) /C 
di sp ( er ) 


Scilab  code  Exa  1.A.7  Additional  Solved  Problem  7 

t=10/1000 
v3=8*exp ( -50*  t ) 
p3  =  v3  * ( -3 . 6) 

P=  [230/1000*4 . 8 ;  ( -3 . 2) *1 . 45 ; p3] 

disp(”Power  absorbed  are  respectively”) 

disp (P) 
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Scilab  code  Exa  1.A.9  Additional  Solved  Problem  9 

Im=250/ 1000 
w  =  100 
r  =  4 

L=50/ 1000 
t=25/1000 


i=Im*sin (w*t) 
Vr=i*r 

Vl=Im*L*w*cos ( w*t ) 


Pr=i*i*r 

Pl=Vl*i 

wl=l/2*L*i*i 

wl_min=0  //min  current  =0 

wl_max = 1 /2* L* Im* Im  //Im  is  the  max  current 

d i s p  ( ”  A 1 1  values  in  Joules”) 
disp(wl_max  ,wl_min  ,  wl  ,P1  ,Pr) 


Scilab  code  Exa  1.1  Problcml 

disp(”vl=A  cos  (w  t+p)”) 

disp  (”  i  1=B  cos  ( wl  t+pl)”) 

disp(”i2=C  sin  (w2  t+p2  )  ”  ) 

d i s  p  ( ” A 1 1  values  in  degrees”) 

w= input ( ”w=” ) 

p= input ( ” p=” ) 

wl  =  input  (”  wl— ’  ) 

pi  =  input (”  pl=” ) 

w2= input (” w2=” ) 

p2= input (” p2=” ) 


if ( w==wl ) 
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lagl=p-pl ; 

disp(lagl,  ”il  lags  vl  by  ”) 

else 

disp(”Lag  undefined  for  il”) 

end 


if  ( w==w2 ) 

lag2=p-p2+90 ; 

disp(lag2,  ”i2  lags  vl  by  ”) 

else 

disp(”Lag  undefined  for  i2”) 

end 


Scilab  code  Exa  1.2  Problem2 

//  X  cos  (w  t+p)+Y  sin  (w  t+q) 


X  =  50 
Y  =  -30 
w  =  120 
p  =  -45 

p  =  p/180*  °/0pi  ; 
q  =  1 60 

q=q/  180*  °/0pi  ; 

/ /  Part  1 

A=X*cos(p)+  Y*sin(q); 
B=-X*sin(p)+  Y*cos(q); 

di sp (B  ,  ” B=” ,A,”A=”) 
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Scilab  code  Exa  1.3  Problcm3 


Gl  =  input (”  Gl=” ) 

G2  =  input (”  G2=” ) 

G3= input (”  G3=” ) 

I s 1 = input (” I s 1 =” ) 

Is2  =  input  (”  Is2=”  ) 

A= [G1+G3 , -G3 ; -G3 , G2+G3] ; 
1= [Isl ; I s2 ]  ; 

V= inv (A) *1 ; 

di sp ( V) 


Scilab  code  Exa  1.5  Problem5 


Rl  =  input (”  Rl=” ) 

R2  =  input (”  R2=” ) 

R3  =  input (”  R3=” ) 

Vs  1  =  input (”  Vsl=” ) 

Vs2  =  input (”  Vs2=” ) 

R= [ R1 +R3 , -R3 ; -R3 , R2+R3] ; 
V  = [Vsl ; -Vs2]  ; 

1= inv (R) *V ; 

di sp ( I ) 


Scilab  code  Exa  1.7  Problcm7 


Rl  =  l 
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R2  =  6 
R3  =  1 2 
R4  =4 
14=3/4 
V6  =6 


V4=3/4*4 
112=3/12 
IR=3/ 4+112 

R= - 12* 1 1 2  +  V6  //KVL 
di sp (R  ,  ”R(01im)=”  ) 

I6=V6/6 

I1=I6+IR 

Vs=l*Il+V6  //KVL 
di sp ( Vs ,  ”  Vs  (V)=” ) 

disp  (Vs*Il  ,  ”  Power  (W)=”  ) 


Scilab  code  Exa  1.9  Problcm9 


i  1  =4 
v3  =  3 
v4  =  8 


i3=v3/3 

i4=v4/4 

i2=(il-i3-i4)/2  //KCL 
v2  =2* i2 


disp (v2 , ” v2=”  ) 
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Scilab  code  Exa  1.10  ProblcmlO 


vl=6 
i2  =2 
13=4 


v2  =2* i2 
v3=2* i3 

v4=-vl+4*i2+v3-v2  / /KVL 

i4=v4/3 

di sp ( i4 ) 


Scilab  code  Exa  1.11  Problemll 

i 1=4-1  //  KCL  at  node  1 
disp(il  i  1=”  ) 


vl  =  i 1  *  1 
vl2=l*l 
v2  =  vl - vl2 
i2=v2 / 1 


i3=l-i2  //KCL  at  node  2 

Vs=v2-l*i3  //KVL 
di sp ( Vs  ,  ”  Vs=”  ) 


Scilab  code  Exa  1.12  Problcml2 
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i34=10-8  //KCL  at  node  4 

v34  =  5* i34 

v23  =40 -10  / /KVL 

vx= v23 

disp (vx , ” vx=”  ) 

ix=4-8  //KCL  at  node  1 

di sp  ( ix  ,  ”  ix=”  ) 

i23=ix+10  //KCL  at  node  2 

R2=vx/ i23 

disp (R2 ,  ”R2=”  ) 

vl4=40+6*ix  //KVL 

Rl=vl4/8 

disp  ( R1  ,  ” Rl=”  ) 


Scilab  code  Exa  1.13  Problcml3 
vs=0.01  //  *cos(1000*t) 
vpi = vs 

i0=-vpi /1000 
vo  =  iO  *  1 000 ; 

printf(”v0  =  %f  *  cos ( 1000* t )”, vo ) 


Scilab  code  Exa  1.15  Problcml5 


t=0 . 1 
R  =  800 

i=50*exp ( - 10*  t ) /1000 
v  =  50*  cos  (25*t ) 
v_i  =  10*t  ~ (2 . 5) 


Pl=i*i*R 

P2=v*v/R 
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9  P3=v_i 
10 

11  disp (P3 , P2 , PI ) 
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Chapter  2 

Fundamentals  of  Resistive 
Circuits 


Scilab  code  Exa  2.1  Probleml 

function  s=series (rl , r2) 
s=rl+r2 
endf unction 

function  p=parallel (rl , r2 ) 
p=rl*r2/ (rl+r2) 
endf unction 

rl=series (12 ,8) 
r2=parallel (20, rl) 
r3  =  series (r2  ,50) 
r4=parallel (30, r3) 
r5=series (10 , r4) 
r6  =  series (r5  ,20) 
Req_ab=parallel (r5 ,40) 
di sp ( Req_ab ) 

r7=40+20+10  //series 
Req_bc=parallel(r4,r7) 
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di sp ( Req_bc ) 


Scilab  code  Exa  2.2  Problem2 

I=(14-4) /C8+5+5+7)  / /KVL 
di sp ( I ) 


Scilab  code  Exa  2.3  Problcm3 

rl =100 
v  =  3/4 

r2=r 1 *v/ ( 1 - v) 


function  p=parallel (rl , r2 ) 
p=rl*r2/(rl+r2) 
endf unction 

////// part  1 

R2_eq  =  parallel (r2  ,  10000) 
k=R2_eq/ (R2_eq+rl) 
change  1  =  (3/4 -k) /( 3/4) *100 
di sp ( change  1 ) 

III  III  Part2 

R2_eq  =  parallel (r2  ,  1000) 
k=R2_eq/ (R2_eq+rl) 
change 2  =  (3/4 -k) /( 3/4) *100 
di sp ( change2 ) 
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Scilab  code  Exa  2.4  Problem4 


function  p=parallel (rl , r2 ) 
p=rl*r2/ (rl+r2) 
endf unction 


v  =  150 

Req  =  1000*parallel (parallel (10 ,5)  ,4) 

il=v/ (2000+0 . 2* 1000+Req) 

vl  =  i 1 *0 . 2*1000 

VReq=Req* i 1 

i2=VReq/ 4000 

disp (vl ) 
disp(i2) 


Scilab  code  Exa  2.5  Problem5 

function  [Rab , Rbc,Rca]=star_to_delta(Ra,Rb,Rc) 
X=Ra*Rb+Rb*Rc+Rc*Ra 
Rab=X/Rc 
Rbc=X/Ra 
Rca=X/Rb 
endf unction 

function  p=parallel (rl , r2 ) 
p=rl*r2/(rl+r2) 
endf unction 

[Rx ,Ry , Rz]  =  star_to_delta (8 ,4 ,2) 

Req_ad=parallel (parallel (4,Rx)+parallel (12 ,Rz) , 
parallel (3 , Ry ) ) 
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di sp ( Req_ad ) 


Scilab  code  Exa  2.6  Problcm6 

//Converting  Voltage  source  with  series  resistance 
to  current  source 

1=1/ (1/4) 

//KCL  at  nodes  1,2,3  using  conductances 

A  =  [11  ,  -3  ,  -4;  -3 ,6  ,  -2;  -4  ,  -2 , 11] 

V= inv (A) *  [5; -7 ; 6] 

di sp ( V) 


Scilab  code  Exa  2.7  Problcm7 

function  p=parallel (rl , r2 ) 
p=rl*r2/ ( r 1 +r2 ) 
endf unction 

Req  =  parallel (15 ,30) 

/ /KCL  at  nodes  1  and  2 

A=[3  ,  — 1;  — 3 ,5] 

V= inv (A) *  [80; -360] 

disp (V  (1)  ) 


Scilab  code  Exa  2.8  Problcm8 
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//Converting  practical  voltage  source  to  current 
source 

1=12/6 

/  /KCLs 

A  = [0 . 7 , -0.5; -0.5 ,0.767] 

V= inv (A) *  [7; 2] 
di sp ( V) 

I12  =  (V (1) -V  (2) ) /2 
disp ( 112) 


Scilab  code  Exa  2.9  Problcm9 

//KCLs 

A=  [2,-1;  -4, 19] 

V= inv (A) *  [10 ; 25] 

Pc  =  l*V  (1) 

Iv= (5-V ( 2 ) ) / 2  +  (5-V  ( 1 ) ) / 5 
Pv=5* Iv 


di sp ( V) 
disp  (Pv  ,  Pc) 


Scilab  code  Exa  2.10  ProblcmlO 

//converting  9v  voltage  source  to  current  source 

1=9/0.45 
1  =  1-7 


//KCLs 
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A  =  [1/0 . 45  +  1/10 , -1/10  ;  0.2-1/10,1/10-0.2] 
V=inv (A) * [ - 18 ; 5] 

disp (V  (2)  ) 


Scilab  code  Exa  2.11  Problemll 

//converting  voltage  source  to  current  source 

1=240/3 

/  /KCLs 

A= [1/3  + 1/6, -1/6; -1/6,  1/6  + 1/1 2  + 1/30] 

V=inv (A) *  [70 ; 5] 

Vx=V (1) 

Vy=V (2) -60 

V6  =  V (2) -V  (1) 

P6=V6*V6/6 

disp (P6) 


Scilab  code  Exa  2.13  Probleml3 
//  KVLs 

A=  [24  ,  -6  ,  -8;  -6 ,20  ,  -5;  -8  ,  -5,25] 
I  =  inv (A) *  [16 ; 0 ; 0] 

15  =  1 (2) -I  (3) 
di sp ( 15 ) 
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Scilab  code  Exa  2.14  Probleml4 


//convert  practical  current  source  to  voltage  source 

V=10*5 

/  /  KVLs 

A  = [35 , -20 ; -20 , 50] 

1  =  inv  (A)  *  [50;  -100] 

v  =  20* ( I (1) -I  (2) ) 
disp (v) 


Scilab  code  Exa  2.15  Probleml5 
11=0.37 

//KVL  for  loop  2 

i2=(-100+20*il) /50 

//KVL  for  loop  1 
R=(50+20*(i2-il) )/ il 
disp (R) 


Scilab  code  Exa  2.16  Problcml6 

//Converting  current  sources  to  voltage  sources 
//  0=[ i  1  ;  i2  ;  vx] 

A=[10  ,  -3 ,4;  -3 ,7 ,0;3  ,  -3  , -1] 
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0  =  inv  (A)  *  [16;  -6  ;  0] 


di sp ( 0 ) 


Scilab  code  Exa  2.17  Probleml7 

//convert  current  sources  to  voltage 

A  = [95 , -15; -15  ,  135] 

I=inv (A) * [20; -70] 

13  =  1 (1) -I  (2) 

160=0.5-1(1) 

V60=I60*60 

P 1 =  0.5*V60  / / O . 5  A  source 

180  =  1  +  1  (2) 

¥80=180*80 
P2=l* V80 

disp ( P2 , PI) 


Scilab  code  Exa  2.18  Probleml8 

//Nodal  Equations 


A=[0. 6,-0. 5 ; 2 , -1.6] 
V= inv (A) *  [4 ; 0] 

V12  =  V (1) -V  (2) 
disp (V12) 
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Scilab  code  Exa  2.19  Problcml9 


//convert  dependent  current  source  to  dependent 
voltage  source 

//Mesh  Equations 

A=[14  ,  -2 ,0;  -2 , 18 ,3;  2  ,  -2  ,  -1] 

0  =  inv (A) *  [  1 00 ; 0  ;  0] 

disp (0  (3)  ) 


Scilab  code  Exa  2.20  Problcm20 


//Nodal  Method 

vl  =  (8/5+16/10) /  (1/5  +  1/  2  +  1/10) 
disp (vl ) 
i5=(8-vl) /5000 
disp(i5) 


//Mesh  Method 

A= [7 , -2; -2 , 12] 

I=inv(A)*[8; -16]/1000 
di sp ( I ) 

vl=2*1000*( I (1) -I  (2) ) 
disp (vl ) 
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Scilab  code  Exa  2.21  Problem21 


function  p=parallel (rl , r2 ) 
p=rl*r2/ (rl+r2) 
endf unction 

//deactivate  voltage  source 

ill =4*6/8 
121=4*2/8 
vl 1 = il 1 *2 
v2 1  =  1  *  1 

//deactivate  current  source 

Req=parallel (2 ,7) 

v22=8*Req/ (2+Req) 
vl2  =  v22*(2/(2  +  5)  ) 

vl=vll+vl2 
v2=v21 +v22 

disp (v2 , vl ) 


Scilab  code  Exa  2.22  Problem22 

function  p=parallel (rl , r2 ) 
p=rl*r2/(rl+r2) 
endf unction 

//Thevenin  Equivalent 
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I=(32-8) / 30 
Voc=32-20*I 
Ro  =  parallel (20 , 10) 
disp (Ro , Voc) 

//Norton  Equivalent 

Isc=32/20+8/10 
disp (Ro , Isc ) 


Scilab  code  Exa  2.23  Problem23 

//convert  current  source  to  voltage 

V=3 . 5*20 

/ /KVL 

1=100/ (10+20+V) 

Voc=100-10*I 

//Finding  Isc 

1=100/10 
Isc  =  ( 1  +  3 . 5) *1 

Ro=Voc/Isc 

disp (Ro , Voc ) 


Scilab  code  Exa  2.24  Problem24 

V  =  32 
1  =  1 
R1  =20 
R2  =8 
R3  =  2 


source 
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Vocl=V/ (R1+R2+R3) *(R2+R3)  /// a,b  open 

Iscl=V/Rl  / / / a , b  short 

Voc2=I*R2/ (R1+R2+R3)*R1 
Isc2=I*R2/ (R2+R3) 

////by  superposition 

Voc=Vocl+Voc2 

Isc=Iscl+Isc2 

disp(Voc) 
disp ( Isc ) 


Scilab  code  Exa  2.25  Problem25 

function  [Ra,Rb,Rc]=delta_to_star (Rab , Rbc , Rea) 
X=Rab+Rbc+Rca 
Ra=Rab*Rca/X 
Rb=Rab*Rbc/X 
Rc=Rca*Rbc/X 
endf unction 

function  p=parallel (rl , r2 ) 
p=rl*r2/(rl+r2) 
endf unction 

[Rl ,R2 , R3]=delta_to_star (4 , 12 ,8) 
Req=10+4+parallel (Rl+15 , R2+16) 

1=12/ Req 
di sp ( I ) 


Scilab  code  Exa  2.26  Problem26 
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//convert  current  source  to  voltage  source 


VI =4*2 
V2  =  2*2 

//Mesh  Analysis 

A=  [7,-3;  -15 , 19] 

1  =  inv (A) *  [6 ; -14] 

Vx  =3*  ( I  (1)  -I  (2)  ) 
disp (Vx) 


Scilab  code  Exa  2.27  Problcm27 

//Short  Circuit  AB  and  convert  Curr  Source  to  V 
source 


V  =  9*5 

Isc=(45-9) / 6 

function  p=parallel (rl , r2 ) 
p=rl*r2/(rl+r2) 
endf unction 

Ro=parallel (5+1 ,3) 
disp (Ro , Isc ) 


Scilab  code  Exa  2.28  Problem28 

//Mesh  Analysis 

A  = [4 , -2;998 ,24.5] 
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1  =  inv (A) *  [1/ 1000 ; 0] 


di sp ( I ) 

P1=I (2) "2*2 . 5*1000 
P2  =  1/1000*1  (1) 
P3=-10~6*I (1) *1 (2) 
P=[P1 , P2 , P3] 
disp (P) 


Scilab  code  Exa  2.29  Problem29 

function  p=parallel (rl , r2 ) 
p=rl*r2/(rl+r2) 
endf unction 

//as  seen  from  ab 

Vthl=30*60/ (30+60) 

Rthl =parallel (30 ,60) 

//as  seen  from  cd 

Vth2=80*40/ (40+60) 

Rth2  =  parallel (40 ,60) 

Vnet  =60  + Vthl -Vth2 
Rnet  =  Rthl  +  Rth2 

disp(Rnet ,Vnet) 


Scilab  code  Exa  2.30  Problcm30 

function  p=parallel (rl , r2 ) 
p=rl*r2/(rl+r2) 
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endf unction 


R=parallel  (60  ,  120) 

//Mesh  Analysis 

A= [6 , -4 ; -4 , 12] 

I=inv (A) * [2.4; -3.6] 
13  =  1 (1) -I  (2) 

160=13*120/ (120+60) 
disp ( 160) 


Scilab  code  Exa  2.31  Problcm31 

//Nodal  Equation 

Vn=198/10 . 5 
I60=Vn/60 
disp ( 160) 


Scilab  code  Exa  2.32  Problcm32 
V8=2*8 
/ /KVL 

Il=(30-16) /10 

//KCL  at  left  node 

12=2-11 


//KVL  middle  mesh 

I4=(10*I1-2*I2) /4 
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Vx=4* 14 


// right  node 

Ix  =  14  - 12 

disp (Vx) 


Scilab  code  Exa  2.33  Problem33 

//////////////////////////////// investigate  the  if 
then  else  construct 


1  =  5 


/ /KVL 

R=(-12.5+15)/I-0.02-0. 035 
disp (R) 

P  =  poly ( [-25 , 12.5  ,0.035]  ,  ”  I”  , ” c o e f f ” ) 
Z=roots (P) 

Ib  =  Z  (2) 
disp(Ib) 

I=(13-12 . 5) /0 . 035 
R= (15-13) /I-O  .  02 
disp (R) 


Scilab  code  Exa  2.35  Problcm35 

R=12* 12/4-25 
disp (R) 
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I  =  sqrt  (1 . 6/1000/10000) 
R=12/I-10000-15000 
disp (R) 


1  =  0 
V  =  12 
di sp ( V) 


Scilab  code  Exa  2.36  Problem36 

function  p=parallel (rl , r2 ) 
p=rl*r2/ ( r 1 +r2 ) 
endf unction 

Req  =  1  +  parallel ( par allel ( 1  ,  1 )  ,1) 
disp (Req) 


Scilab  code  Exa  2.37  Problcm37 

function  p=parallel (rl , r2 ) 
p=rl*r2/(rl+r2) 
endf unction 

Rl  =  parallel (parallel (40 , 60)  ,80+40) 

Req=20+parallel (10+R1 ,30) 
disp (Req) 
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Scilab  code  Exa  2.38  Problem38 


function  p=parallel (rl , r2 ) 
p=rl*r2/ (rl+r2) 
endf unction 

Rl =parallel (5 ,40) +15 
Req=parallel (Rl , 15) +10 
i =6/ ( 1000*  Req) 

ix=i*Rl / (Rl +15) 
di sp ( ix ) 

p  =  i*i*10000  ////// wrongly  done  in  the  book  as  ix*ix 

*10000 

disp (P) 


Scilab  code  Exa  2.39  Problem39 

function  p=parallel (rl , r2 ) 
p=rl*r2/(rl+r2) 
endf unction 

Rl  =  parallel (1  ,3) 

//convert  current  source  to  voltage 

V=20/1000*2000 
I = (40 -8) / (2+0. 75) 


Pr  =  8*  I 
di sp ( Pr ) 


Scilab  code  Exa  2.40  Problcm40 
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function  p=parallel (rl , r2 ) 
p=rl*r2/ ( r 1 +r2 ) 
endf unction 

11=12/1000 
Ix=Il*30/ (30+60) 
di sp ( lx ) 

12=20/1000 
V2  =  12  * ( 10  +  5) 

Rl  =  10  +  parallel  (30 ,60) 

I1=V2/R1 

Ix=Il*30/ (30+60) 

disp(Ix)  //Wrongly  printed  in  the  book  as  10 

Ix=6/1000 

II = (30+60) /30* lx 

V2=I1*R1 

I2=V2 / 15 

disp(I2) 

I  s  =45 

11=45* (10+5) / ( 10+5+R1 ) 

Ix=Il*30/ (30+60) 

disp(Ix)  ///Answer  is  wrong  in  the  book 
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Chapter  3 

Fundamentals  of  Reactive 
Circuits 


Scilab  code  Exa  3.1  Probleml 

function  C=seriesC (Cl , C2) 
C=C1*C2/(C1+C2) 
endf unction 

Ceql=10+seriesC  (10,10) 

Ceq2=Ceql 

Ceq=seriesC ( seriesC (Ceql  ,Ceq2)  ,10) 
disp (Ceq) 


Scilab  code  Exa  3.2  Problem2 

function  p=parallel (rl , r2 ) 
p=rl*r2/ (rl+r2) 
endf unction 
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Leq  =  l  +  parallel (2 ,2) +parallel (parallel (3 ,3)  ,3) 


disp (Leq) 


Scilab  code  Exa  3.3  Problcm3 

function  C=seriesC (Cl , C2) 
C=C1*C2/(C1+C2) 
endf unction 

V  =  100 

Cp=l+2 

Ceq=seriesC (Cp  ,6) 
q=Ceq* V 

ql=q*(l/(l+2) ) 
q2=q-ql 


disp  (q2  ,  ql  ,  q) 


E6=q*q/2/6/1000 
El=ql*ql/2/l/1000 
E2=q2*q2/2/2/1000 
Enet  =  E1  +  E2  +  E6 

disp(Enet ,E6 ,E2 , El ) 


Scilab  code  Exa  3.4  Problem4 

L  =  poly([3,-4,l]  , ”L”  , ” c  o  e  f f  ”  ) 
disp (roots (L) ) 
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Scilab  code  Exa  3.5  Problem5 


function  C=seriesC (Cl , C2) 
C=C1*C2/(C1+C2) 
endf unction 

Ceq  =  l+seriesC (6 ,3) 

disp(Ceq) 


Scilab  code  Exa  3.6  Problcm6 

function  p=parallel (rl , r2 ) 
p=rl*r2/ (rl+r2) 
endf unction 

Lbc  =  parallel  (1/10  ,  1/20) 

Leq=l/50+Lbc 

disp (Leq) 

Vbc=Lbc*150 
di sp ( Vbc ) 


Scilab  code  Exa  3.7  Problcm7 

function  C=seriesC (Cl , C2) 
C=C1*C2/(C1+C2) 
endf unction 
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Ceq= ser i e sC ( ser ie sC  ( 10 , 20)  ,40) 
disp(Ceq) 


q=Ceq*280 
disp (q) 

VI =q/10 
V2=q/20 
V3=q/40 

disp(V3 ,V2 , VI) 


Scilab  code  Exa  3.8  Problem8 

Ceq=10+20+40 
Q=200* [10 ,20 ,40] 

disp (Q) 

Qnet=Q (1) +Q (2) +Q (3) 
di sp ( Qnet ) 


Scilab  code  Exa  3.9  Problcm9 

//deactivate  all  independent 

function  p=parallel (rl , r2 ) 
p=rl*r2/(rl+r2) 
endf unction 


L  =  2 

Req=parallel (6 , 12) +4 

T1=L/Req 

disp ( T 1 ) 


sources 


46 


11 

12 

13 

14 

15 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

1 

2 

3 

4 

5 

6 


C  =  l/2 

Req  =  2  +  parallel (6 , 12) +2 

T2=Req*C 

disp (T2) 


Scilab  code  Exa  3.11  Problcmll 


L  =  2 


VL_0minus=0  //steady  state 

disp(VL_0minus) 

i_0minus=l 
i_0plus=i_0minus 
VL_0plus=12-4*i_0plus 
di sp ( VL_0plus ) 

di_by_dt_0plus=VL_0plus/L 

disp(di_by_dt_Oplus) 


Scilab  code  Exa  3.12  Problcml2 

Vc_0m=100 
Vc_0p=¥c_0m 
T= (400+100) *2E-6 
i_0p=100/500 
P_0p  =  i_0p  "2*400 
di sp ( P_0p ) 
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Scilab  code  Exa  3.13  Problcml3 


Vc_0m=10 
Vc_0p=Vc_0m 
disp ( Vc_0p ) 

/ /KVL 

il_0p=(10-Vc_0p) /20 
i2_0p=Vc_0p/20 

//KCL 

iC_0p=il_0p-i2_0p 
disp ( iC_0p ) 

iC_inf=0  //capacitor  is  open  circuit 

disp ( iC_inf ) 

VC_inf =10*20/ (20+20) 
disp ( VC_inf ) 


Scilab  code  Exa  3.16  Probleml6 

L  =  0 . 25 
R  =  250 
V  =  10 
T  =  L/R 

t=0 . 5E-3 

i=V/R*(l-exp(-t/T)) 
di sp ( i ) 

t=2E-3 

i=V/R*((exp((lE-3-t)/T)-exp(-t/T))) 
di sp ( i ) 
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Scilab  code  Exa  3.18  Probleml8 


R=5000 
C=1 . OE-6 
Ie=l/1000 
Vcf =Ie*R 
T  =  R*  C 


t=10/1000 

VclO=Vcf*(l-exp(-t/T)) 

disp(VclO) 

t=20/1000 

Vc20=Vcf*(l-exp (-t/T) ) -VclO 
disp ( Vc20 ) 


Scilab  code  Exa  3.21  Problem21 

//Replace  curr  source  by  voltage 

Ics=l/1000 
R=100* 1000 
V= Ics  *R 


Vc_0m=25/ (25+100) *V 
disp ( Vc_0m) 

Vc_0p=Vc_0m 
disp ( Vc_0p ) 

Vc.inf =(V-10) *25/ ( 100+25) 
disp ( Vc_inf ) 


Scilab  code  Exa  3.22  Problem22 


source 
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L=20/ 1000 
VR2_0p=10 
VR2_inf  =0 
VR2_ 1 =5 

t=l/1000 

T=-t / log ( (VR2_1 -VR2_inf ) / ( VR2_0p - VR2_inf ) ) 
R2=L/T 

Rl=l/(2/VR2_0p -1/R2) 
disp (R2 , R1 ) 


Scilab  code  Exa  3.24  Problem24 
C  =  l/4 

Vc_0m=12/ (12+6) *12 
Vc_0p=8 
disp ( Vc_0p ) 
t  =  0 

V_0p  =  6*  cos (t) 

//at  t=0+ 

Vth=12/ ( 12+6) *6 
Rth=6*12/ (6+12) 
ic_0p=(Vth-8) /Rth 
disp(ic_0p) 
d_by_dt_Vc_0p=ic_0p/C 
disp(d_by_dt_Vc_Op) 


Scilab  code  Exa  3.26  Problem26 

//at  t  =0— 
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iL=4 
v  =  0 

disp(v,iL) 

//at  t=0+ 

iL=4 
v= -4*20 
disp(v,iL) 


Scilab  code  Exa  3.27  Problem27 

function  p=parallel (rl , r2 ) 
p=rl*r2/(rl+r2) 
endf unction 

L=25/ 1000 

//at  t  =0— 

Rl  =  parallel (150 ,75) 

iL_0m=30/2 

disp ( iL_0m) 

R2  =  parallel (150 , 50) 
ix_0m=R2/ (R2+75) *30 
disp(ix_0m) 


//at  t  =0+ 

iL_0p= iL_0m 

T  =  L/  (75  +  50) 
iL_inf  =0 

//at  t  =0.2  ms 

t  =0 . 2E-3 
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iL=iL_inf+(iL_Op-iL_inf)*exp(-t/T) 

ix=-iL 

disp(ix,iL) 


Scilab  code  Exa  3.29  Problem29 

R  =  250 
C=25E-3 


//at  t=0— 

v_0m=200*0 . 2 

//at  t=0+ 

v_0p=v_0m 


T  =  R*C 
v_inf  =0 

//at  t  =0.2  ms 

t=0 . 2E-3 


t=T*log((v_0p-v_inf)/(20-v_inf)) 
disp (t) 


Scilab  code  Exa  3.32  Problcm32 


L  =  0 . 2 
R=50+30 

iL_0m= 100/50 
iL_0p= iL_0m 
iL_inf  =0 
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T  =  L/R 


t=10/1000 

iL_10=iL_0p*exp(-t/T) 
disp ( iL_10) 

tl=-T*log(0.5*iL_0p/(iL_0p)) 
disp  (tl ) 


Scilab  code  Exa  3.33  Problem33 

R  =  5 
L  =  2 
V  =  10 
T  =  L/R 

iL_m=V/R 
disp(iL_m) 
iL_0p= iL_m 
iL_inf  =2  +  V/R 

//function  I=i  L  ( t ) 

//  I=iL_inf  +  (iL_0p  —  iL_inf)  *exp(— t  /T) 
/  /  endfunction 

//  disp  ( iL  ( 0 . 5 )  ) 

//disp  (iL  (1.5)  ) 


t  =0 . 5 

I=iL_inf+(iL_Op-iL_inf)*exp(-t/T) 
di sp ( I ) 

t  =  l  .  5 

I=iL_inf+(iL_Op-iL_inf)*exp(-t/T) 
di sp ( I ) 
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Scilab  code  Exa  3.36  Problcm36 

L=4/ 1000 
R=1000 
V  =  9 


iL_0m=V/ (2*R) 
iL_0p= iL_0m 
iL_inf  =0 
isw_inf  =V/R 
T  =  L/R 

t=5E-6 

iL_5=iL_0p*exp(-t/T) 

disp(iL_5) 

isw_5=isw_inf*(l-exp(-t/T)) 

disp(isw_5) 


Scilab  code  Exa  3.37.38  Problcm37  and  38 


R  =  20 
L  =  2 

C=l/25 

//at  t=0— 

iR  =  -4 
iL=4 
iC  =  0 
vR=-R*4 
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vC=-vR 
vL  =  0 


disp(iR,  iL  ,  iC  , vR , vL  ,vC) 

//at  t=0+ 

iL  =  iL 
vC  =  vC 

iR=5-4  //KCL  at  node  1 
vR  =  R* iR 

iC=4+iR  //KCL  at  node  2 
vL  =  vC  +  vR  //KVL  inner  mesh 

disp(iR,  iL  ,  iC  , vR , vL  ,vC) 


//at  t  =0+  derivatives  are 

DiL=vL/L 

DvC=iC/C 

DiR  =  -DiL  //Differentiating  KCL  at  node  1 

DiC  =  DiR  //Differentiating  KCL  at  node  2 

DvR  =  R*  DiR 

DvL  =  DvR  +  DvC  //Differentiating  KVL  inner  mesh 

disp(DiR ,DiL ,DiC ,DvR ,  DvL ,DvC) 
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Chapter  4 

Steady  State  Analysis  For 
Sinusoidal  Excitation 


Scilab  code  Exa  4.1  Probleml 

1  lav  =  1/  (2*  °/„pi  )  *  int  egrat e  (  ’  20*  sin  ( wt )  ’  ,  ’  wt  ’  ,  %pi  /3  ,  %pi  ) 

2  disp ( lav) 

3 

4  Irms  =  sqrt(l/(2*°/„pi)*integrate(,(20*sin(wt))~2’,  ’wt’, 

°/„pi  /3  ,  °/„pi  )  ) 

5  disp ( Irms ) 


Scilab  code  Exa  4.2  Problem2 

1  Iav_by_Im  =  l/°/0pi*integrate(  ’sin(wt)  ’  ,  ’wt’  ,0,°/0pi) 

2  disp ( Iav_by_Im) 

3 

4  //  final  answer  is  wrong  in  the  book 
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Scilab  code  Exa  4.4  Problem4 


/////////////////////Part  a 

L  =  0 . 05 
R  =  20 
w=l . 2E3 
Xl=w*L 

V=complex  (100 ,200) 

I  =  complex  (5,0) 

Xc=-imult (R+imult (Xl)-V/I) 
C=l/w/Xc 

disp (C) 


/////////////////////Part  b 

w  =  200 

V=complex  (100 ,0) 

Xl=w*L 

Xc=-imult (R+imult (Xl)-V/I) 
C=l/w/Xc 

disp (C) 


Scilab  code  Exa  4.5  Problem5 

L  =  0 . 0255 
Rl=6 
V  =  240 
f  =50 

w  =  2*  °/0pi  *f 
Va=sqrt (V~2/10) 

Vb=3* Va 
Xl=w*L 

Thetal=atan(Xl/Rl) 
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Theta2  =  7„pi /2-Thetal  //V a  and  Vb  are  in  quadrature 


I = Vb / sqrt (R1 ~2+Xl ~2) 

R=(Va/I)/ sqrt (l+tan(Theta2) "2) 
disp (R) 

Xc=R*tan(Theta2) 

C=l/w/Xc 
disp (C) 


Scilab  code  Exa  4.6  Problem6 

Zl  =  complex  (10 , 15) 

Z2=complex (6 , -8) 

1  =  complex (15 ,0) 

I1=I*Z2/ (Z1+Z2) 

I2=I*Z1/(Z1+Z2) 

phase=[atan(imag(Il)/real (II) ) ; atan(imag(I2)/real (12 

))] 

disp(phase*180/%pi) 

V=I1*Z1 

disp (180/%pi*atan(imag(V)/real (V) ) ) 


Scilab  code  Exa  4.7  Problem7 

function  [r , theta] =cart_to_polar (z) 
x=real (z) 
y= imag (z) 
r=norm (z) 
theta=atan(y/x) 
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endf unction 


function  Zeq=parallel ( Z 1 , Z2) 
Zeq=Zl*Z2/ (Z1+Z2) 
endf unction 

w  =  400 

R  =  1 0 

L=25E-3 

C=250E-6 

Xl=w*L*%i 

Xc  =  l/  ( w  *  C  *  °/0  i  ) 


//////////// Part  a 

Zinl =R+Xc 

[rl ,thetal]=cart_to_polar(Zinl) 
disp(thetal*180/%pi ,rl) 


//////////// Part  b 

// Zin2=10+  p  ar  a  1 1  e  1  (Xc  ,  XI )  /// impedence  is  infinite 

...  thus  an  error 

//  [r2  ,theta2]=cart_to_polar  (Zin2) 

disp(90,”  inf”) 


////////////Part  c 

Zin3=R+parallel (Xc , Xl+10) 

[r3 , theta3] =cart_to_polar (Zin3) 
disp(theta3  *180/70pi  ,r3) 


Scilab  code  Exa  4.8  Problcm8 

function  [r , theta] =cart_to_polar (z) 
x=real (z) 
y= imag (z) 
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r=norm (z) 
theta=atan(y/x) 
endf unction 

function  [x,y]=polar_to_cart (zpolar) 
r=real (zpolar ) 
theta=iraag  (zpolar)  /180*°/0pi 
x  =  r *  cos (theta ) 
y=r*sin (theta) 
endf unction 

function  Zeq=parallel (Z1 , Z2) 
Zeq=Zl*Z2/ (Z1+Z2) 
endf unction 

R  =  200 
L  =  0 . 5 
C=50E-6 
w  =  200 

Xl=w*L*%i 
Xc  =  l/  (w*C*°/0i) 


////////////////////////////Part  a 

Ir  =  complex(0.02*cos(30*°/„pi/180)  ,0.02*sin(30*°/„pi/180) 

) 

V= Ir *R 
I1=V/X1 
Ic=V/Xc 
I=Ir+Il+Ic 
di sp ( I ) 


////////////////////////////Part  b 

[lx  ,  Iy] =polar_to_cart ( complex (2 , -40)  ) 
I=complex (lx , Iy) 

Zin=parallel (parallel (R , Xc) ,X1) 

V=Zin* I 
Ir=V/R 
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disp(V,Ir) 


Scilab  code  Exa  4.9  Problcm9 

p  =  poly(  [0,2, 0,-8]  ,”w”  ,”  coeff”) 
w=roots (p) 
disp(w(l)  ) 


Scilab  code  Exa  4.10  ProblcmlO 

function  Zeq=parallel ( Z 1 , Z2) 
Zeq=Zl*Z2/ (Z1+Z2) 
endf unction 

Zeq  =  parallel (5 , 3- imult  (6) ) 

I  =  complex  (10  ,0) 

V=Zeq* I 

pf=cos (at an ( imag (V) / real (V) ) ) 

P=norm(V)*norm(I)*pf 

disp (pf  ,  P) 


Scilab  code  Exa  4.11  Problcmll 

PI =8000 
pf 1=0 . 8 
V 1 =430 

Il=Pl/Vl/pf 1 
Prl=Pl-Il ~2*1 
Ql=Pl*tan(acos (pfl) ) 
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Qr 1  =  Q1 - 1 1 ~  2* 1 . 2 
Srl=sqrt(Prl~2+Qrl~2) 

Vl  =  Sr 1 / 1 1 
PI =15000 

Ql=Pl*tan(acos (0.8) ) 
Pr2=Pl -Pr 1 
Qr2=Ql -Qr 1 

Sr2  =  sqrt (Pr2  ~2+Qr2  ~2) 

1 2  =  Sr2  /  VI 
P2=Pr2+I2~2*0 .7 
Q2  =  Qr2  + 12 "2*0 . 9 

pf  =  cos (atan(Q2/P2)  ) 

S2  =  sqrt (P2  ~2+Q2  ~2) 
V2=S2/I2 

disp(V2 ,pf , P2 ,Q2  ,Q1 ,  VI) 


Scilab  code  Exa  4.12  Problcml2 

V  =  23 1 
f  =50 

w  =  2*  °/0pi  *f 
Sl  =  10 
pf =0 . 8 
Pl=Sl*pf 

Ql=Sl*sin(acos (pf) ) 
Qc=Ql-tan(acos (0.95) )*P1 

C=Qc*1000/V"2/w 

Isl=10*1000/V 

Il=Isl 
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Is2  =  norm(Pl  +  0/0i*(Ql-Qc)  )  *1000/V 


disp(Is2 , Isl ,C,Qc) 


Scilab  code  Exa  4.13  Probleml3 

function  [x,y]=polar_to_cart (r, theta) 
theta  =  theta/180*°/0pi 
x  =  r *  cos (theta ) 
y=r*sin (theta) 

endf unction 

//converting  v  source  to  curr  source 

V  =  15 

R  =  5 

1 1  =  V/R 


[Ix,Iy]=polar_to_cart (2,-30) 

I2=complex (lx , Iy) 

/ /  nodal  analysis 

A=  [0.7-0.  l*°/0i  ,  -0.5;  -0.5 ,0.5  +  0.5*%i] 

V= inv (A) *  [II  ;  12] 

[rvl  ,thetavl]=polar(V(l)) 

[rv2 ,thetav2]=polar(V(2)) 

disp(real(thetav 2*180  / °/0pi)  ,  rv2,real(thetavl*180/°/„pi) 
,  rvl ) 


Scilab  code  Exa  4.14  Problcml4 
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w  =  600 
R  =  1 0 
L=10E-3 
C=0 . 5E-3 
Xl=w*L*%i 

Xc  =  l/(w*C*°/0i)  //Xc  value  in  the  book  is  wrong 

V 1  =40 

V2  =  complex (0  ,  -  30) 

/ / mesh  analysis 

A= [R+Xl ,  -R ; -R , R+Xc] 

1  =  inv  (A)  *  [VI  ;  V2] 

Ir  =  I (1) -I  (2) 

Vr =  R* Ir 


disp(Vr)  //Final  answer  is  different 


Scilab  code  Exa  4.15  Probleml5 

Zl=complex (2 , -1) 

Z2=complex (0 , -5) 

Z3=complex (2,1) 

11  =  1 

I2=complex (0 , -2) 

//deactivate  curr  source  2 

V11=I1*(Z2+Z3)/(Z1+Z2+Z3)*Z1 
V21=I1*Z1/ (Z1+Z2+Z3)*Z3 

//deactivate  curr  source  1 

V12=I2*Z3*Z1/ (Z1+Z2+Z3) 
V22=I2*(Z1+Z2) / (Z1+Z2+Z3) *Z3 

V1=V11+V12 
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V2=V21+V22 


disp (V2 , VI) 


Scilab  code  Exa  4.16  Problcml6 

function  [Ra,Rb,Rc]=delta_to_star (Rab , Rbc , Rea) 
X=Rab+Rbc+Rca 
Ra=Rab*Rca/X 
Rb=Rab*Rbc/X 
Rc=Rca*Rbc /X 
endf unction 

function  Zeq=parallel ( Z 1 , Z2) 

Zeq=Zl*Z2/ (Z1+Z2) 
endf unction 


V  =  20 
w  =  2 
R  =  1 
C  =  l/2 

Xc  =  l/  ( w  *  C  *  °/0  i  ) 

Zab=-%i 
Zbc  =2 
Zca  =  2 

[Za,Zb ,Zc]=delta_to_star(Zab ,Zbc , Zca) 
Zth=Za+parallel (R+Zb , Xc+Zc ) 

Vdc=V* (Zc+Xc) /(R+Zb+Zc+Xc) 
di sp ( Zth , Vdc ) 

Isc=Vdc/Zth 
disp (Zth , Isc ) 
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Scilab  code  Exa  4.17  Probleml7 


function  Zeq=parallel ( Z 1 , Z2) 
Zeq=Zl*Z2/ (Z1+Z2) 
endf unction 


V  =  20 

Zl  =  complex (5 , 10) 

Z2=complex (3 , -4) 

Vth=V*Z2/ (Z1+Z2) 

Zth=parallel ( Z 1 , Z2) 

//by  maximum  power  transfer  theorem 

Zl=conj (Zth) 

P=norra(Vth/ (Zth+Zl) ) ~  2*  real (Zl) 
disp(P,Zl) 

Rl=sqrt (real (Zth) ~2+(4+imag(Zth) ) ~  2) 
disp (R1 ) 


Scilab  code  Exa  4.18  Problcml8 

function  [x,y]=polar_to_cart (r, theta) 
theta  =  theta/180*°/0pi 
x  =  r *  cos (theta ) 
y=r*sin (theta) 
endf unction 


Rl  =  l 
R2  =  6 
L  =  2 
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C  =  l/3 


[Vx,Vy]=polar_to_cart  (4,30) 

[lx , Iy]=polar_to_cart (0.8 , -10) 

//deactivate  current  source 

V=complex(Vx ,Vy) 
w  =  3 

Xl=w*L*%i 
X  c  =  1  /  (w*C*°/0i) 

Ill=V/(Rl+R2+Xl+Xc) 
Pavl=norm(Ill) ~2*R1 

//deactivate  voltage  source 

I=complex (lx , Iy) 
w  =  5 

Xl=w*L*%i 
Xc  =  l/  ( w  *  C  *  °/0  i  ) 


112=1* (R2+X1) /(R2+Xl+Rl+Xc) 
Pav2=norm(I12) ~2*R1 

Pav=Pavl+Pav2 

disp(Pav) 


Scilab  code  Exa  4.21  Problem21 


1=35 
VR  =  25 
Vrl =40 
VRr 1 =50 
V  c  =45 
C=50E-6 
Xc  =  ¥c / 1 
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w=l/(Xc*C) 


theta=  acos ( (VR~2+VRrl ~2-Vrl ~2) / (2* VR*VRrl ) ) 

x=VRrl*cos (theta) -25 

y=VRrl*sin(theta) 

r  =  x/I 
L=y/ ( I *w) 

Vappl=sqrt ((VR+x) ~2+y ~2) 

R=VR/I 

di sp  (L  ,  r  ,  R) 


Scilab  code  Exa  4.22  Problcm22 


V  =  12 
w=1000 
R  =  4 

L=5E-3 
C=400E-6 
Xl  =  w*L*°/>i 
Xc  =  l/  ( w  *  C  *  °/0  i  ) 


/ / mesh  analysis 

A= [R+Xl , -XI ; 2 -XI , Xl+Xc] 

I  =  inv (A) *  [12 ; 0] 

disp(I(D)  //answer  is  wrong  in  the  book 


Scilab  code  Exa  4.23  Problem23 

function  [x,y]=polar_to_cart (r, theta) 
theta  =  theta/180*°/0pi 
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x  =  r *  cos (theta ) 
y=r*sin (theta) 
endf unction 

[I lx , Ily]=polar_to_cart (2  ,  10) 
[I2x , I2y]=polar_to_cart  (3 , 120) 

Il=complex (Ilx , Ily) 

I2=complex (I2x , I2y) 

w  =  4 

R  =  2 

L  =  3 

C  =  l/4 

Xl=w*L*%i 

Xc  =  l/  ( w  *  C  *  °/0  i  ) 


//deactivate  source  2 

il=(R+Xl) /(R+Xl+Xc) *11 
i2=l/ (R+Xl+Xc) *12 
i=il+i2 

disp(i)  //answer  is  wrong  in  the  book 


Scilab  code  Exa  4.24  Problem24 
V  =  230 

w  =  2*  °/0pi  *50 
R  =  8 
L  =  0  .  1 
C=160E-6 

Xl=w*L*%i 
disp (XI ) 

Xc  =  l/  ( w  *  C  *  °/0  i  ) 
disp(Xc) 
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Z=norm(R+Xl+Xc) 
disp (Z) 

I  =  V/Z 
di sp ( I ) 

pf=cos (atan(norm(Xl+Xc) /R) ) 
disp (pf ) 

Vcoil=I*norm(R+Xl) 

Vc=norra ( I *Xc ) 
disp (Vc , Vcoil ) 


Scilab  code  Exa  4.25  Problem25 

function  Zeq=parallel ( Z 1 , Z2) 
Zeq=Zl*Z2/ (Z1+Z2) 
endf unction 

V  =  100 

Zl=complex (0 , -5) 

Z2  =  complex  (5,0) 

Z3  =  complex (0,5) 

I  =  [0 ,0,0] 

Z=Z1 +parallel (Z2 , Z3) 

I ( 1 ) =V/Z 

I  (2)  =  Z3  /  (  Z2  +  Z3  )  *  I  (1) 

I (3) =Z2/(Z2  +  Z3)*I  (1) 

di sp ( I ) 


Scilab  code  Exa  4.26  Problem26 
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function  Zeq=parallel ( Z 1 , Z2) 
Zeq=Zl*Z2/ (Z1+Z2) 
endf unction 

//converting  curr  source  to  v 

1=0.4 
R1  =  15 
V=I*R1 
w  =  400 


R2  =5 
R3  =  80 
L=25E-3 
Xl=w*L*%i 


Zeq=Rl+R2+parallel (R3 ,X1) 
I=V/Zeq 

11=1* R3/ (R3+X1) 

Ix=I*Xl/ (R3+X1) 

disp  (II  ,  lx) 


Scilab  code  Exa  4.27  Problem27 


V  =  4 
w  =  2 


R  =  4 
L  =  2 

Xl=w*L*%i 


I 1 = V/ (Xl+R) 
Xc=-V/ imag ( II ) 
C= 1 / (w*Xc) 
disp (C) 


source 
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Scilab  code  Exa  4.28  Problem28 


function  [x,y]=polar_to_cart (r, theta) 
theta  =  theta/180*°/0pi 
x=r*cos (theta) 
y=r*sin (theta) 
endf unction 


1=4 
w  =  500 

[Ir , Ic]=polar_to_cart (3 ,40) 
Ix=complex (Ir , Ic) 

C=lE-3 
L=10E-3 
Xl=w*L*%i 
Xc  =  l/  ( w  *  C  *  °/0  i  ) 


Il=I-Ix 
Ic  =  lx 
Vc  =  Ic  *Xc 
V1=I1*X1 
Vx=Vl - Vc 

disp(¥x)  //answer  wrong  in  the  book 


Scilab  code  Exa  4.29  Problem29 


w  =  10 
R  =  2 
L  =  0 . 3 

Ir=10*(cos  (°/„pi/4)+°/0i*sin(%pi/4)  ) 
Vr  =  Ir  *R 
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Vc=20*(cos  (°/„pi/4)  -°/„i  *  sin  ( %pi  /  4)  ) 

V=Vr+Vc 
Ic  =  Ir 

C=Ic/ (w*Vc*%i) 

di sp  ( real  ( C)  )  ///////Wrong  in  book 


Scilab  code  Exa  4.30  Problcm30 

R=2000 

C=5E-6 

w=1000 


Xc  =  l/  (  w  *  C  *  °/0  i  ) 

Y= 1 / R+ 1 / Xc 
Z  =  1  /  Y 

Req=real (Z) 

Ceq=-1/ (imag(Z)*w) 

disp (Ceq*lE6 , Req)  ////// Answer  wrong  in  the  book 


Scilab  code  Exa  4.31  Problem31 

1  =  10*  (cos  (°/0pi/180*37)  -°/„i*sin(°/0pi/180*37)  ) 
V  =  6 

C=250E-6 

w=1000 

Xc  =  l/  (  w  *  C  *  °/0  i  ) 


Ic=V/Xc 

disp(Ic) 

Ix= imag ( I - Ic ) 
X=abs (V/Ix) 
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Ir=real ( I - Ic ) 
R=V/Ir 

disp (R , X) 


Scilab  code  Exa  4.32  Problem32 

Vl=30  +  10*°/„i 
V2  =30 
w=1000 
L  =  1 

C 1 = 1 E -6 
C2  =  1 E -6 
R1 =1000 
R2  =  1 000 


Xl=%i*w*L 
Xcl  =  °/„i/(w*Cl) 

Xc2  =  °/0i  /  ( w*C2 ) 

//////// mesh  equations 

A  =  [  1  ,  % i  ;  °/„i  ,l-2*'/.i] 

I  =  inv  (A)  *[10*°/„i  ;  30] 

Icl  =  I (1) -I  (2) 

disp(Icl) 


Scilab  code  Exa  4.33  Problcm33 

function  [x,y]=polar_to_cart (r, theta) 
theta  =  theta/180*°/0pi 
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x  =  r *  cos (theta ) 
y=r*sin (theta) 
endf unction 

[Ir , Ic]=polar_to_cart (3 ,30) 
I=complex(Ir , Ic) 

Z  c  =  -  5  *  °/0  i 
Zrl  =  10  +  5*°/0i 

Ic=Zrl/ (Zc+Zrl)*I 
Vth=Ic*Zc 

Zth=parallel (Zc , Zrl) 
In=Vth/Zth 

disp(Rth , In ,Vth) 


Scilab  code  Exa  4.34  Problem34 

function  [x,y]=polar_to_cart (r, theta) 
theta  =  theta/180*°/0pi 
x  =  r *  cos (theta ) 
y=r*sin (theta) 
endf unction 

[Ir , Ic] =polar_to_cart (2 , 10) 

Il=complex (Ir , Ic) 

[Ir , Ic]=polar_to_cart (3 ,90  +  30) 

I2  =  complex (Ir  ,  Ic) 


w  =  4 
R  =  2 
L  =  0 . 5 
C  =  l/4 
Xl=w*L*%i 
Xc  =  l/  ( w  *  C  *  °/0  i  ) 
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//deactivate  voltage  source 

il=Il*(R+Xl)/ (R+Xl+Xc) 

//deactivate  curr  source 

i2=I2/(R+Xl+Xc) 

i=il+i2 
di sp ( i ) 


Scilab  code  Exa  4.35  Problcm35 

function  Zeq=parallel ( Z 1 , Z2) 
Zeq=Zl*Z2/ (Z1+Z2) 
endf unction 


1  =  10 
w=2000 
R  =  200 
L  =  0 . 125 

Xl=w*L*%i 


Zeq=parallel (R , XI ) 

V  =  I  *  Zeq 
I1=V/X1 
Ir=V/R 

t=l/1000 

I I  =  sqrt (2) *  real ( I  *  exp ( °/„i  *  w*t ) ) 

III  =  sqrt(2)*real(Il*exp(°/„i*w*t)) 
Ir 1  =  sqrt (2) *  real ( Ir *  exp (%i*w*t) ) 

VI  =  sqrt  (2)  *real  (V*exp  (°/0i*w*t)  ) 


Ps  =  - VI  *  1 1 
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Pr=-Vl * Ir 1 
P1=-V1*I11 

disp  (PI  ,  Pr  ,  Ps  ) 


Scilab  code  Exa  4.36  Problcm36 

function  [x,y]=polar_to_cart (r, theta) 
theta  =  theta/180*°/0pi 
x=r*cos (theta) 
y=r*sin (theta) 
endf unction 

[Vr ,Vc]=polar_to_cart (120 ,45) 
V=complex(Vr ,Vc) 

[Zr,Zc]=polar_to_cart (16.3 ,24.5) 
Z=complex(Zr ,Zc) 


w  =  50 


I  =  V/Z 

P  =  V*  con j  ( I ) 

Pavg=real (P) 

pf=real ((V/I)/norm(V/I) ) 
di sp (pf , Pavg ) 


Scilab  code  Exa  4.37  Problcm37 

function  Zeq=parallel (Z1 , Z2) 
Zeq=Zl*Z2/ (Z1+Z2) 
endf unction 
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function  [x,y]=polar_to_cart (r, theta) 
theta  =  theta/180*°/0pi 
x=r*cos (theta) 
y=r*sin (theta) 
endf unction 

[Ir , Ic]=polar_to_cart  (20 ,60) 
I=complex(Ir , Ic) 

w=5000 
R=3000 
L  =  1 

C=0 . 25E-6 
Xl=w*L*%i 
Xc  =  l/  ( w  *  C  *  °/0  i  ) 


Z=parallel (R+Xl , Xc) 
V  =  I*Z 
di sp ( V) 


Scilab  code  Exa  4.38  Problem38 

function  Zeq=parallel (Z1 , Z2) 
Zeq=Zl*Z2/ (Z1+Z2) 
endf unction 


V  =  5 


w  =  5 
R  =  2 

L 1 =0 . 8 
L2  =  1 
C  =  l/5 

Xll=w*Ll*°/„i 
X12  =  w*L2*°/„i 
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Xc  =  l/  (w*C*°/0i) 


Vth=R/ (R+X12)*V 
Zth=Xc+parallel (X12 ,  R) 
di sp ( Zth , Vth ) 

I=Vth/ (Zth+Xll) 
V1=I*X11 
S  =  V1 *  con j  (I) /2 
di sp ( S ) 


Scilab  code  Exa  4.39  Problcm39 
Pm=600 

Vm=sqrt (3*Pm) 
f =15 . 36E6 
w  =  2*  °/0pi  *f 
t  =20 . 3E-3 

theta  =  yopi/2-modulo  ( w*t  ,  2 * °/„p i  ) 
disp(theta*180/%pi) 


Scilab  code  Exa  4.40  Problem40 

V  =  240 
w  =  100*°/0pi 
R  =  6 


Vr=120 

I=Vr/R 

t  =  (205/I)~2  ////t=rA2+XP2 
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r = ( (240/ I ) “2-t -R*R) /2/R  ///this  part  solved 

wrong  in  the  book 
Xl  =  sqrt (t-r*r) 

Z=sqrt (t) 
disp (r) 
disp (XI ) 
disp (Z) 

Pl_choke=I*I*r 
disp (Pl_choke ) 
pf=Pl_choke/205/20 
disp (pf ) 


Scilab  code  Exa  4.41  Problem41 

w  =  2 
L  =  1 
C  =  0 . 5 
Xl=w*L*%i 
Xc  =  l/  (  w  *  C  *  °/0  i  ) 


V  =  1 

I=V/(l  +  */.i) 

Y  =  I  /  V 

R=l/ real (Y) 
di  sp (R , Y) 


Scilab  code  Exa  4.42  Problem42 

w  =  400 
R  =  5 

L=25E-3 
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C  =  1 . 25E-3 
Xl=w*L*%i 
Xc  =  l/  (w*C*°/0i) 


Z=R+X1 +Xc 
Y  =  1  /  Z 

C=- imag ( Y) /w 
Yn=real ( Y) 
Rn=l/Y 

disp (C) 


Scilab  code  Exa  4.43  Problcm43 

function  Zeq=parallel ( Z 1 , Z2) 
Zeq=Zl*Z2/ (Z1+Z2) 
endf unction 

w  =  800 

R1  =  10 

R2  =20 

L  =  0 . 05 

C=0 . 25E-3 

Xl=w*L*%i 

Xc  =  l/  ( w  *  C  *  °/0  i  ) 


Zeql=R2+parallel (Xc , Rl+Xl) 
Zeq2=parallel (parallel (R1 ,R2)+Xc ,  XI ) 

disp ( Zeq2 , Zeql ) 


Scilab  code  Exa  4.44  Problem44 
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pf =0 . 707 
theta=  acos (pf ) 

P  =  200 

Q=tan (theta) *P 
pf 2=0 . 85 

theta_n=acos (pf2) 
Qn=Q*tan(theta_n) 
Qc=Q-Qn 


w  =  3 14 
V= 2000 

C  =  Qc/ (V*V*w) 
disp (C) 


Scilab  code  Exa  4.45  Problem45 


1=20 
w=2000 
R  =  200 
L  =  0 . 25 
Xl=w*L*%i 


Ir=I*Xl/(Xl+R) 

11=1 - Ir 
V1=X1*I1 

t=lE-3 

ir  =  sqrt(2)*real(Ir*exp(°/„i*w*t)) 
il  =  sqrt  (2)  *  real  (Il*exp(°/0i*w*t)  ) 
vl  =  sqrt  (2)  *  real  (Vl*exp(°/0i*w*t)  ) 
is  =  sqrt  (2)  *real  (I*exp  (°/„i*w*t)  ) 
vs  =  vl 

Pr=ir*ir*R 
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Pl=vl* il 

Ps=is*ir*R 

Pr=ir*vl 


disp  (Ps  ,  PI  ,  Pr ) 


Scilab  code  Exa  4.46  Problem46 

Vl  =  l 
w  =  2 
R  =  1 
r  =  l/2 
L  =  0 . 5 
C  =  1 

Xl=w*L*%i 
Xc  =  l/  ( w  *  C  *  °/0  i  ) 


I1=V1/X1 
VR= II  *R 
Vs=¥l+VR 

Vr  =  Vs 
Ir=Vr/r 
Ic=Vs/Xc 
Is=Ir+Ic+Il 

PR=- VR* VR/R 
Pr=- Vr * Vr /r 
Ps  =  Vs  *  conj ( Is ) 

disp(real(Ps) , Pr , PR) 


Scilab  code  Exa  4.47  Problem47 
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V  =  1 0  *  %  i 
w=1000 
R  =  2 

L=2E-3 
C=500E-6 
Xl=w*L*%i 
X  c  =  1  /  (w*C*°/0i) 


Zeq=Xl+Xc+parallel (R , XI) 
disp (Zeq) 

I  =  V/Z 

Ir=I*Xl/(R+Xl) 

112=1 - Ir 

V1=I*X1 

Vc=I*Xc 


disp  (I , Ir , 112 , VI , Vc) 


Scilab  code  Exa  4.48  Problem48 

function  [x,y]=polar_to_cart (r, theta) 
theta  =  theta/180*°/0pi 
x=r*cos (theta) 
y=r*sin (theta) 
endf unction 

[Vr,Vc]=polar_to_cart  (100,120) 
Vl=complex (Vr , Vc ) 


1  =  10 
w=5000 
R=10000 
L  =  2 . 4E-3 
C=0 . 05E-3 
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[Vr ,Vc]=polar_to_cart (96 ,60) 
Vl=complex(Vr,Vc) 


Xl=w*L*%i 
X  c  =  1  /  (w*C*°/0i) 


I1=V1/X1 

Ic=I-Il 

V2=(R+Xc) *  I  c  +  V 1 
V3=V2 - VI 
Vl=Ic*Xc 

disp (VI , V2 , V3) 


Scilab  code  Exa  4.49  Problcm49 

function  Zeq=parallel (Z1 , Z2) 
Zeq=Zl*Z2/ (Z1+Z2) 
endf unction 


////////////short  AB 

Vab  =0 
Isc=50/4 


////////////for  Zn 

Zn  =  parallel  (4 ,  8*  °/0i ) 

disp  ( Isc  ,  Zn) 


Scilab  code  Exa  4.50  Problcm50 
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function  Zeq=parallel ( Z 1 , Z2) 
Zeq=Zl*Z2/ (Z1+Z2) 
endf unction 


V  =  1 2  *  %  i 


Vth  =  4-12*°/0i  /  (4-12*y„i+6  +  9*°/0i)  *V 
Zth  =  parallel  (4 -12* % i  ,6  +  9*70i) 

I  =  Vth/ ( Zth  +  6  +  12* %i ) 

S  =  V*  con j  ( I ) 
di sp ( S ) 

Zl  =  conj  (Zth) 

I=Vth/ ( Zth+Zl ) 

S  =  V*  con j  ( I ) 
di sp ( Zth , S ) 


Scilab  code  Exa  4.51  Problem51 

function  Zeq=parallel ( Z 1 , Z2) 
Zeq=Zl*Z2/ (Z1+Z2) 
endf unction 


1=5 


Zth  =  parallel  (40 ,  3  0  *  7,  i  ) 
Zl=conj (Zth) 

Vth=Zth* I 

Z=Zth+Zl 

Il=Vth/Z 

Praax=abs (Il*Il*real (Zl) ) 


di  sp  ( Zth  ,  Praax ) 
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Scilab  code  Exa  4.52  Problem52 


V=4*  sqrt (2) 
w  =  2 
R=4 
L  =  2 

Xl=w*L*%i 


1 1  =  V/ ( R  +  Xl ) 

Ic  =  - imag ( II )  * °/„i 
C  =  Ic/(V*°/0i*w) 


disp (C) 
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Chapter  5 

Frequency  Response 


Scilab  code  Exa  5.1  Probleml 
Q0  =  200 

//a 

R=  1 

C=2E-6 

L=C/ ( QO/R) ~ 2 
disp (L) 


L=2E - 15 
C  =  1 . 2E-9 
R=Q0*  sqrt (L/C) 
disp (R) 

R=1 18500 
L=120E-12 
C  =  L* (Q0/R) ~ 2 
disp (C) 


Scilab  code  Exa  5.2  Problem2 
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R=1000 

C=49E-6 

L=13E-3 


w0=l/ sqrt (L*C) 
Q0  =  w0  *R*C 

di  sp  (  Q0  ,  wO  ) 


Scilab  code  Exa  5.3  Problcm3 
R  =  5 

C=0 . 08E-6 
L=5E-3 


w0=l/ sqrt (L*C) 
f  0  =  w0 /2/ Zpi 
Q0  =  w0  *L/R 
bandwidth=wO/QO 
wb=bandwidth 
w=[w0-l/2*wb ,w0+l/2*wb] 

disp(w(2)  ,  w  (1)  ) 


Scilab  code  Exa  5.4  Problem4 


L=40E-6 
R=4 . 02 
f  0  =  800  E3 


wO  =  2* %pi *f  0 
C  =  1 / ( w0*w0*L) 
Q0=w0*L/R 
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fb=fO/QO 


Xl  =  wO*L*°/,i 

Xc=l/(wO*C*%i) 

Zin=R+Xl+Xc 

di  sp  ( Zin  ,  f  b  ,  C ) 


Scilab  code  Exa  5.5  Problcm5 


G=5E-6 
L=2E-3 
I =2E -3 
w=5000 

C=l/ (w*w*L) 
Vmax=I/ G 
11  =  1 

I=Il*w*L*G 
di sp ( I , Vmax , C ) 


Scilab  code  Exa  5.6  Problcm6 

w0=1000 
w2  =  1 050 
wb  =  2* (w2-w0) 

Q0=w0/wb 
di sp ( Q0 ) 
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Scilab  code  Exa  5.7  Problcm7 


wO  =  200  E3 
wb=5E3 
L  =  2 . 5E-3 
Q  =65 


Q0=w0/wb 
r=wO*L/ Q 

disp (r) 


Scilab  code  Exa  5.8  Problem8 


C=12E-9 
L=4E-3 
R  =  5 


w0=l/ sqrt (L*C) 

Vc  =  1 . 5 
I=wO*C*Vc 
Zi  =  R 
Vi=R*I 


disp(Vi ,w0) 


Scilab  code  Exa  5.9  Problem9 

wO  =  2 . 5  E6 
Zin=60E3 
Q0  =80 
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R=60E3 
C=QO/(wO*R) 
L=l/ (C*wO*wO) 

disp(C,L,R) 


Scilab  code  Exa  5.10  ProblcmlO 


C  =  2 . 5E-6 
R  =  8 

wb=400 
L=R/ wb 

w0=l/ sqrt (L*C) 

Q0  =  w0  *L/R 
wl =w0 - wb /2 
w2=w0+wb/2 

disp(w2 ,wl ,Q0 ,w0 ,L) 


Scilab  code  Exa  5.11  Problemll 

w  =  1 E6 
wb=1000 
V=0 . 05 
I=5E-3 

R=V/I 
Q0=w0/wb 
L  =  R*  q0 / wO 
C=l/ ( w0*w0*L) 

Vl=wO*L*I 
Vc  =  ¥  1 
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wl =wO - wb /2 
w2=w0+wb/2 

disp(w2 ,  wl ,Vc ,V1 ,C,L,Q0 ,  R) 


Scilab  code  Exa  5.12  Probleml2 

R=  1  E6 
L  =  1 

C=lE-6 

I=10E-6 


w0=l/ sqrt (L*C) 
V  =  I  *R 


wb=l/ (R*C) 

Q0=w0/wb 
wl =w0 - wb /2 
w2=w0+wb/2 

disp(w2 ,  wl ,  V  ,  Q0 ,w0) 


Scilab  code  Exa  5.13  Probleml3 
f 0=1E6 

Craax=500E-12 
C=450E-12 
wO  =2*  %pi *f 

L=l/(wO*wO*Cmax) 

w=l/ sqrt (L*C) 
f  =w/  (2*  °/„pi  ) 
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wb  =  2*2*°/„pi*(f-f0) 
r=wb*L 

Q0=2*%pi*f*L/r 


disp  ( Q0  ,  L  ,  r) 

////////calculation  mistakes  in  book 


Scilab  code  Exa  5.14  Probleml4 

R=10E3 

L=50E-3 

C=100E-6 

w0=l/ sqrt (L*C) 

Q0  =  w0  *R*C 

alpha= wO /2/ Q0 

zeta=alpha/wO 

wd  =  w0  *  sqrt (l-zeta~2) 

disp(wd ,zeta , alpha , Q0 ,w0) 

//////////calculation  mistakes  in  the  book 


Scilab  code  Exa  5.15  Problcml5 

w0=1000 
wd=997 
Yin=l . 2E-3 

zeta=sqrt (l-(wd/w0) ~2) 

alpha=zeta*wO 

R=l/Yin 
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C  =  1 / (2* alpha*R) 
L=l/ (w0*w0*C) 

disp(C,L,R) 


Scilab  code  Exa  5.17  Problcml7 


1  =  1 

w  =  400 

R  =  500 

L=l/40 

C=250E-6 

Xl=w*L*%i 

Xc  =  l/  ( w  *  C  *  °/0  i  ) 


w0=l/ sqrt (L*C) 
////////at  resonance 

Ir  =  I 
V  =  R*  I 
I1=V/X1 
Ic=V/Xc 
Icir=abs ( II ) 

Ic+Il=0 

Icl=Ic+Il 

disp(Icl , Ic , II , Ir) 


Scilab  code  Exa  5.18  Probleml8 

V  =  100 
R  =  1 0 
L=2E-3 
C=200E-6 
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w0=l/ sqrt (L*C) 

Xl  =  wO*L*°/,i 
Xc=l/(wO*C*%i) 

I  =  V/R 

V1=I*X1 

Vc=I*Xc 

Vlc=Vl+Vc 

disp (Vic , Vc ,V1 ,1 ,w0) 


Scilab  code  Exa  5.19  Probleml9 


R  =  1 

L=10E-6 
C=10E-12 
V  =  10 


w0=l/ sqrt (L*C) 

Xl  =  wO*L*°/0i 
Xc=l/(wO*C*%i) 

Q0=w0*L/R 
Z  =  R*  Q0  *  Q0 
Iin=V/Z 

Ic=V/Xc 

disp(Ic , Iin ,Z,Q0 ,w0) 


Scilab  code  Exa  5.20  Problcm20 
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V  =  200 
R  =  2 

L  =  0 . 02 


f  =25 

wO  =2*  %pi *f 
C  =  1 / ( wO*wO*L) 
I  =  V/2 

Vc=l/ (C*wO) 
di sp ( I , Vc ) 


f  =50 

wO  =2*  %pi *f 
C  =  1 / ( wO*wO*L) 
I  =  V/2 

Vc=l/ (C*wO) 
di sp ( I , Vc ) 

f  =  100 

wO  =2*  %pi *  f 
C  =  1 / ( wO*wO*L) 
I  =  V/2 

Vc=l/ (C*wO) 
di sp ( I , Vc ) 


Scilab  code  Exa  5.21  Problem21 


R  =  1 0 
L  =  0  .  1 
C=150E-6 

w=sqrt ((1+R*R*C/L)/ (L*C) ) 
f =w/ (2*  %pi ) 

di sp ( f  ) 
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Req=R/(l-w*w*L*C+(w*R*C) ~  2) 
disp (Req) 


Scilab  code  Exa  5.22  Problem22 


R  =  35 
Q0  =50 
f 1=540E3 
f 2=1610E3 
wl =2*  %pi *f 1 
w2=2*%pi*f 2 


L=l/(wl*(Q0/R) ) 

Craax=(Q0/R) ~2*L 
Cmin  =  l/ ( L  *  w2  ~  2 ) 

disp(Cmin*1000000 , Cm ax *1000000 ,L*1000000) 


Scilab  code  Exa  5.28  Problem28 

function  [x,y]=polar_to_cart (r, theta) 
theta  =  theta/180*°/0pi 
x  =  r *  cos (theta ) 
y=r*sin (theta) 
endf unction 

[Vr,Vc]=polar_to_cart (4,30) 
V=complex(Vr ,¥c) 

[Ir , Ic]=polar_to_cart (0.8 , -10) 
I=complex(Ir , Ic) 
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wv  =  3 
wi  =5 
R1  =  1 
L  =  2 
R6  =6 
C  =  l/3 

//////////deactivating  curr  source 

Xl  =  wv*L*°/0i 
Xc  =  l/  (wv*C*°/0i) 

Ill=¥/(R1+R6+Xl+Xc) 

Pavl=norm(Ill) ~2*  R1 

//////////deactivating  voltage  source 

Xl  =  wi*L*°/0i 
Xc=l/ (wi*C*%i) 

112=1* (R6+X1) /(R6+Xl+Rl+Xc) 
Pav2=norm(I12) ~2*R1 

Pav=Pavl+Pav2 

disp(Pav) 


Scilab  code  Exa  5.29  Problem29 


L=10E-6 
R  =  1 

C=10E-9 
V  =  10 


Zraax=L/R/ C 
I=V/ Zmax 
w0=l/ sqrt (L*C) 
w=0 . 9* wO 
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11  Y  =  R*C/L  +  0/0i*(w*C-l/w/L) 

12  I  =  norra  (Y)  *V 

13 

14  disp(I) 
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Chapter  6 

Three  Phase  Circuits 


Scilab  code  Exa  6.1  Probleml 

function  [x,y]=polar_to_cart (r, theta) 
theta  =  theta/180*°/0pi 
x=r*cos (theta) 
y=r*sin (theta) 
endf unction 

omega  =  exp(°/0i*120/180*%pi) 

Va=200 

Vb  =  200/ omega 
Vc  =200* omega 

[Zr , Zc] =polar_to_cart  (100  ,  60) 
Z=complex(Zr ,Zc) 

Ia= Va/Z 
Ib=Vb/Z 
Ic=¥c/Z 

disp  (Va  ,  Vb  ,  Vc  ,  la  ,  lb  ,  Ic) 
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Scilab  code  Exa  6.2  Problem2 


VI =400 
Vph=Vl 
pf =0 . 8 

Pph=1500/3 
Iph=Pph/Vph/pf 
11  =  sqrt (3) * Iph 
theta=  acos (0.8) 

Iph  =  Iph*exp(°/0i*theta) 
Zph=Vph/ Iph 

disp (Iph , II , Zph) 


Scilab  code  Exa  6.3  Problcm3 

VI =400 
pf =0 . 8 

Pph  =  1200/3 
Vph=Vl /sqrt (3) 
Iph=Pph/Vph/pf 
11= Iph 

theta=acos (0.8) 

Zph  =  Vph/Iph*exp(70i*theta) 

disp(Il , Iph , Zph) 
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Scilab  code  Exa  6.4  Problem4 


function  [x,y]=polar_to_cart (r, theta) 
theta  =  theta/180*°/0pi 
x=r*cos (theta) 
y=r*sin (theta) 
endf unction 

VI =400 

Vp=Vl/sqrt (3) 
disp (Vp) 

[Ir , Ic]=polar_to_cart (50 , -30) 
Il=complex (Ir , Ic) 

Ip  =  Il 
Zy=Vp/Il 
di sp ( Zy ) 

P=sqrt (3) *Vl*Il*cos (30/ 180*  %pi ) 

Q  =  sqrt  (3)  *Vl*Il*sin  (30/ 180*  °/„pi ) 
disp (P , Q) 


Scilab  code  Exa  6.5  Problcm5 

VI =400 
Z  =  16  +  °/0i*12 

Vp  =  Vl 
Ip=Vp/Z 

pf=cos (at an ( imag (Z) / real (Z) ) ) 
11  =  Ip  *  sqrt  (3) 

P=sqrt (3) *Vl*Il*pf 
Q  =  sqrt (3) *Vl*Il*sin(acos (pf)  ) 
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S  =  P  +  %  i  *  Q 


disp  ( Ip  ,  II  ,  pf  ,  P  ,  Q  ,  S) 


Scilab  code  Exa  6.6  Problcm6 
VI =400 

Zload  =  60  +  °/„i*15 

Vs=Vl/sqrt (3) 

Z  =  0 . 3  +  %i *  1  +  Zload 

Il=Vs/Z 

Vload=Il*Zload*sqrt (3) 
Pload=3*Il*Il*real (Zload) 
Qload=3*Il*Il*imag (Zload) 

Sload  =  Pload  +  °/0i  *  Qload 

//////////// 1 1  =  1 i n  e 1 o  s  s 

Pll =3* II *11* real (Z- Zload) 
Qll=3*Il*Il*imag(Z-Zload) 

Sll  =  Pll  +  °/.i*qil 

Ssource=Sload+Sll 

disp (II , Vload ,Sload ,S11 , Ssource) 


Scilab  code  Exa  6.7  Problcm7 

Igl=200000/ sqrt (3) /11000/0.75 
PG1 =200000 

QGl=200000*tan(acos  (0.75) ) 

SGl  =  PGl+°/„i*QGl 
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Pll=3*Igl*Igl*l 

Qll=3*Igl*Igl*2.2 

sn=Pii+,/.i*qn 

S1G1 =SG1 -Sll 
V1=S1G1 / sqrt (3) /Igl 


PI =400000 

Ql=400000*tan(acos (0.8) ) 
P1G2=P1 -real (S1G1 ) 
QlG2=Ql-imag(SlGl) 

SlG2  =  PlG2  +  °/0i*qiG2 

IG2  =  S1G2/ sqrt (3)  /VI 

P112=3*IG2*IG2*0 . 6 
qil2=3*IG2*IG2*l . 2 

PG2=P1G2+P112 
qG2=qiG2+qii2 
SG2  =  PG2  +  °/0i*qG2 
VG2  =  SG2 / sqrt (3) /IG2 


di sp ( norm ( SG2 ) ,norm(VG2)) 


Scilab  code  Exa  6.9  Problcm9 

function  [x,y]=polar_to_cart (r, theta) 
theta  =  theta/180*°/0pi 
x=r*cos (theta) 
y=r*sin (theta) 
endf unction 

VI =400 
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[Ir , Ic]=polar_to_cart (20 , -30) 

Iab  =  complex(Ir  ,  Ic) 

omega  =  exp(°/0i*120/180*%pi) 
Ibc=Iab/ omega 
Ica=Iab*omega 

IAa=Iab-Ica 

P ab  =  Vl*norm (I ab ) *  cos (30/180*  °/„pi ) 
Ptotal =3*  Pab 

Zp= VI / lab 
Rp=real ( Zp ) 

disp ( IAa) 
disp (Ptotal ) 
disp (Rp) 


Scilab  code  Exa  6.10  ProblcmlO 

omega  =  exp(°/oi*120/180*%pi) 

Zp  =  1 7 . 3  2  +  %  i  *  1 0 
Vab  =400 
Iab=Vab/Zp 
Ica=Iab/ omega 
IAa=Iab-Ica 
IBb= IAa/ omega 
ICc  =  IAa*  omega 
disp(ICc , IBb , IAa) 

Pab=(norm(Iab) "2) *real (Zp) 
Ptotal =3*  Pab 
disp (Ptotal ) 
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Itotal=IAa+IBb+ICc 
disp ( Itotal ) 


Scilab  code  Exa  6.11  Problemll 

omega  =  exp(°/oi*120/180*%pi) 

function  [x,y]=polar_to_cart (r, theta) 
theta  =  theta/180*°/0pi 
x=r*cos (theta) 
y=r*sin (theta) 
endf unction 

[Zr,Zc]=polar_to_cart (5,60) 
Z=complex(Zr ,Zc) 

[Vr ,Vc]=polar_to_cart (25 ,30) 
Van=coraplex(Vr ,Vc) 

I an = Van /Z 
Ibn= Ian/ omega 
Icn=Ian*omega 

Vcn=Van*  omega 
Vac=Van-Vcn 

disp(Ibn , Icn ,Vac) 


Scilab  code  Exa  6.12  Probleml2 
V  =  400 

w  =  2*50*  °/„pi 
P=25000 
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pf 1=0.7 

thetal=acos (pfl) 

Ill=P/(  sqrt  (3)  *V*pf  1)  *exp  ( -°/„i *thetal ) 
Ipl=Ill/sqrt (3) 


pf 2=0 . 85 

theta2=acos (pf2) 

I12  =  P/(sqrt(3)*V*pf2)*exp(-°/0i*theta2) 

Ip2=I12/ sqrt (3) 

Ic=Ip2-Ipl  ////calculation 

mistake  in  the  book  at  this  step 

C=real (Ic/(¥*w*%i) ) 
disp (C) 


Scilab  code  Exa  6.13  Probleml3 

omega  =  exp(°/0i*120/180*%pi) 

¥rn=400/sqrt (3) 

¥yn=Vrn/ omega 
¥bn=¥rn* omega 

PI =4000 
P2  =8000 
P3  =  1 2000 

Ir=conj (Pl/¥rn) 

Iy=conj (P2/¥yn) 

Ib=conj (P3/¥bn) 

In=Ir+Iy+Ib 

disp(norm(Ir) ,norm(Iy) ,norm(Ib) ,norm(In)) 
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Scilab  code  Exa  6.14  Probleml4 


function  [x,y]=polar_to_cart (r , theta) 
theta  =  theta/180*°/0pi 
x=r*cos (theta) 
y=r*sin (theta) 
endf unction 


Vbc  =  -  400*  %i 
Rl  =  l 


[Sr , Sc]=polar_to_cart (10000 , acos (0.8) ) 
S=complex (Sr , Sc) 


Vl=Vbc 

Il=norm(S/(sqrt(3)*Vl))  ////////calculation 

mistake  in  book  here 
P  =  3* II *  II *R1 
disp (P) 

Van=norm(Vbc/ sqrt (3) ) 

Ia=Il*exp(-%i*acos (0.8) ) 
omega  =  exp(°/0i*120/180*%pi) 

Ib= la/ omega 

disp(Van,Ia,Ib) 


Scilab  code  Exa  6.15  Probleml5 
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function  [x,y]=polar_to_cart (r, theta) 
theta  =  theta/180*°/0pi 
x=r*cos (theta) 
y=r*sin (theta) 
endf unction 

[Vr ,Vc]=polar_to_cart  (231  ,60) 
Van=coraplex(Vr ,Vc) 

Sp  =  ( 2 . 5-*/.i*l  .  2)  *1000 

omega  =  exp(0/0i*120/180*%pi) 

Vbn=Van/ omega 
Vcn  =  Van*  omega 

Vbc=Vbn-Vcn 

Ibc=conj (Sp/Vbc) 

Ica=Ibc*  omega 
Iab= Ibc / omega 

IaA=Ica-Iab 

disp ( IaA) 


Scilab  code  Exa  6.16  Probleml6 

omega  =  exp(7.i*120/180*7.pi) 

Vab  =400 
Vca=400*  omega 
Vbc  =400/ omega 

//Mesh  Method 

A  =  [80  +  100*%i  ,  100*7.1  ;  100*7.i  ,50*'/.i] 
I=inv(A)*[-Vca;Vbc] 
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Ia= I  ( 1 ) 
Ib  =  I  (2) 


Ic=- ( Ia+Ib) 
Van=80* la 
Vbn  =  -50*0/0i*Ib 
Vcn  =  100*°/0i*Ic 


disp(Ia,Ib,Ic,Van,  Vbn  ,  Vcn ) 


Scilab  code  Exa  6.17  Probleml7 

omega  =  exp(°/0i*120/180*%pi) 

Vab  =400 

Vca  =  400*  omega 

¥bc  =400/ omega 

I  ab  =  V  ab  /  ( 1 0  0  *  °/0  i  ) 

Ibc  =  Vbc/(-50*°/„i) 

Ica=Vca/80 

Ia=Iab-Ica  /// c alaculat i on  mistake  in  the  book 

here 

Ib= Ibc - lab 
Ic=Ica-Ibc 

di  sp ( la , lb , Ic ) 


Scilab  code  Exa  6.18  Probleml8 


V  =  400 
P  =  900 
pf =0 . 8 
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Pph=P/3 
Vph=V/ sqrt (3) 
Iph=Pph/Vph/pf 
Zph=Vph/ Iph 
theta=  acos (pf ) 

Zph  =  Zph*exp(-°/„i*theta) 
di sp ( Zph ) 


Scilab  code  Exa  6.19  Probleml9 
V  =  400 

Z  =  complex (40 ,30) 

Iph=¥/Z 

11= sqrt (3) *norm(Iph) 

Ptotal=sqrt (3) *Il*V*cos (at an ( imag (Z) / real (Z) ) ) 
disp (Iph , II , Ptotal) 


Scilab  code  Exa  6.20  Problem20 

Igl =15000/ sqrt (3) /800/0 . 8 
PG1 =15000 

QGl  =  15000*tan(acos  (0.8) ) 
SGl  =  PGl  +  °/0i*QGl 

Pll =3* Igl * Igl *  1 .2 
Qll=3*Igl*Igl*l  .8 

sn=Pii+°/.i*qii 

S1G1 =SG1 -Sll 
V1=S1G1 / sqrt (3) /Igl 
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PI =30000 

Ql=30000*tan(acos (0.8) ) 
P1G2=P1 -real (S1G1 ) 
QlG2=Ql-imag(SlGl) 

SlG2  =  PlG2  +  °/0i*QlG2 

IG2=S1G2/ sqrt (3) /VI 

P112=3*IG2*IG2*0 . 8 
Q112=3*IG2*IG2*1 . 2 

PG2=P1G2+P112 
QG2=QlG2+qil2 
SG2  =  PG2  +  °/0i*qG2 
VG2  =  SG2 / sqrt (3) /IG2 


di sp ( norm ( SG2 ) ,norm(VG2)) 


Scilab  code  Exa  6.21  Problem21 

omega  =  exp(°/„i*120/180*°/opi) 

w=2*50*  °/„pi 

Vry  =4 1 5 

Vyb=Vry/ omega 

Vbr =Vry *  omega 


C=40E-6 
Zry =100 

Zyb  =  complex (20 ,60) 
Zbr =1/ ( %i *w*C) 

1 1  =  Vry / Zry 
I2=Vyb/Zyb 
13= Vbr / Zbr 
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15 

16  Ir  =  I  1  - 13 

17  Iy  =  I2  - 1 1 

18  Ib  =  I3  - 12 

19 

20  di sp ( norm ( Ir ) ,norm(Iy) , norm(Ib)) 


114 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

1 

2 

3 


Chapter  7 

Magnetic  Circuits 


Scilab  code  Exa  7.1  Probleml 

uo  =  (4* °/0pi)  *lE-7 

i  1  =5 
i2  =  2 . 5 
r  =  0 . 4 

H  =  il/(2*°/„pi*r) 

F=uo*H*i2  ///////// attractive 

ur  =  8000 
Firon=ur*F 

di sp (H , F , Fir on ) 


Scilab  code  Exa  7.2  Problem2 

1=0 . 15 
11=50 
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H 1 = i 1 / ( 2 *  %p i *  ( 0  . 1+1) ) 
i2  =  -Hl*(2*0/0pi*0. 1) 


disp(i2) 


Scilab  code  Exa  7.3  Problcm3 

Ha  =  4/  (2*°/0pi*0.2) 

Hb  =  Ha 

H=sqrt (Ha~2+Hb~2) 

theta=(°/„pi  +  atan(-Hb/Ha)  )  *180/°/0pi 
disp (H , theta) 


Scilab  code  Exa  7.4  Problem4 


uo  =  (4*°/„pi)  *lE-7 

Bg= 1 ■ 2 
N  =  400 
ur  =4000 

lc= (2* (20-4+16-4) -0 . 2) /100 

Ac=4/100*4/100 

Rc=lc/(ur*uo*Ac) 

lg  =  0 . 2/100 
Rg= lg / (uo*Ac) 

R=Rc+Rg 

f lux=Bg*Ac 
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i=f lux*R/N 


di sp ( i ) 


Scilab  code  Exa  7.5  Problem5 

function  Zeq=parallel ( Z 1 , Z2) 
Zeq=Zl*Z2/ (Z1+Z2) 
endf unction 


uo  =  (4*7,pi)  *lE-7 
f luxgl=0 . 8E-3 

lgl =0 . 02/100 
lg2  =  0 . 02/100 
lg3  =  0 . 025/100 
Agl=2/100*l/100 
Ag2=l/100*l/100 
Ag3=l/100*l/100 

Rgl=lgl / (uo*Agl ) 
Rg2=lg2/(uo*Ag2) 

Rg3=lg3/ (uo*Ag3) 
Req=Rgl+parallel (Rg2 ,  Rg3) 

MMF  =fluxgl*Req 

disp (MMF) 


Scilab  code  Exa  7.6  Problcm6 


uo  =  (4*°/0pi)*lE-7 
ur  =4000 
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fluxc=0.01 
N  =  500 


LI = (2* (20+4) +25+4) /100 
Lc=(25+4) /100 
Lr  =  L1 -0 . 02/100 
Lg=0 . 02/100 
A=4/ 100*4/ 100 

R1=L1/ (ur*uo*A) 

Rc=Lc/ (ur*uo*A) 
Rr=Lr/(ur*uo*A) 

Rg  =  Lg / ( uo  *  A ) 

Fab  =  f luxe  *Rc 

fluxr=Fab/(Rr+Rg) 

fluxl=fluxc+fluxr 

F=fluxl*Rl+Fab 
I  =  F/N 

di sp ( I ) 


Scilab  code  Exa  7.7  Problcm7 


uo  =  (4*yopi)*lE-7 


11=50/100 
12=20/100 
13=50/100 
lbc  =0 . 025/100 


Al=25E-4 
A2  =  1 2 . 5E-4 
A3=25E-4 
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f luxg=0 . 75E-3 


B=f luxg/ A1 
Fbc  =  B/uo  * lbc 

Hcd=200 

Hab=Hcd 

Fabcd=Hab*ll 

Fad=Fbc+Fabcd 

Had=Fab/12 

Bad=l . 04 

f luxad=Bad*A2 

fluxdea=fluxad+fluxg 
Bdea=fluxdea/A3 
Hdea=500 
Fdea=Hdea* 13 

F=Fdea+Fad 

disp (F) 


Scilab  code  Exa  7.8  Problem8 

uo  =  (4*7,pi)  *lE-7 
ur  =  3000 

Ac=10/ 100*10/ 1000 
Be  =  1 . 4 

Hc=Bc/ (uo*ur ) 
lc=150/ 100 
Ftotal=Hc*lc 

N2  =  800 
12=2 
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F2=N2* 12 
Fl=Ftotal -F2 

N 1 =400 

I1=F1/N1  ////out  of  terminal  a 

disp(Il) 


Be  =  1 .4 

fluxc=Bc*Ac 

Rc=lc/(Ac*uo*ur) 

Wf  =  l/  2*Rc*fluxc*  f luxe 
disp (Wf ) 


L1=N1*N1/Rc 
L2=N2*N2/Rc 
M=  sqrt (L1*L2) 
di sp (M) 


Scilab  code  Exa  7.9  Problcm9 


uo  =  (4*°/„pi)  +1E-7 

ur  =  6000 

A=5/ 2*2/10000 

Lr  =  (  °/0pi  *  (20  +  25)  /2-0.  1)  /100 
Lg=0 . 1/100 
Rr=Lr/(uo*ur*A) 

Rg=Lg/ (uo*A) 

Rtotal=Rr+Rg 

F=2*500 
f lux=F/Rtotal 
Bg=f lux/A 
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disp (Bg) 

L  =  f lux  *500/2 
disp (L) 

Wfr=l/2*Rr*flux*flux 
Wfg=l/2*Rg*flux*flux 
disp (Wf g , Wf r ) 

iraax=2 

f lux_max=f lux 

E  =  4 . 44*314/(2*  °/0pi)*500*flux_inax 

f  lux_max2  =  100/ (  sqrt  (2)*4. 44*314/(2*  0/0pi)*500) 

Fmax2=flux_max2*Rtotal 

imax2=Fraax2/500 

di sp ( imax2 ) 


Scilab  code  Exa  7.10  ProblcmlO 

Ac=12/ 10000 
Fmin=160*9 . 81 
B=sqrt (Fmin*2*uo/Ac) 

H=2800 
L=75/ 100 
F  =  H*L 

Lg=0 . 1/1000 
A=24/ 10000 
Rg  =  Lg/ ( uo  *  A ) 
f luxg=B*A 
Fg=f luxg*Rg 

Ftotal=F+Fg 

irain=Ftotal/1000 

disp(imin) 
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Scilab  code  Exa  7.11  Problcmll 

Wml=l/2* (1200 -750) *0.012 
disp (Wml ) 

f lux_g=0 . 012/750*1200 

Wm2  =  1 /2* ( f lux_g -0.012) *1200 

disp (Wm2) 

Ra=1200/0 .012 
Rb=1200/f lux_g 
Fo=-1/2*(0.012~2) *0.75*10~6 
Fc=-l/2*( f lux_g "2) *0.75*10~6 
disp (Fc , Fo) 


Scilab  code  Exa  7.12  Probleml2 
Vl  =  10 

//Mesh  Method 

A  =  [1  +  10*  %i  ,  - 8 0  *  % i  ;  8 0  *  % i  ,  -500 -800*  %i  ] 
1  =  inv  (A)  *  [VI  ;  0] 

V2  =500* I  (2) 

disp (VI , V2  ,  I) 


Scilab  code  Exa  7.14  Probleml4 
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uo=(4*%pi) *lE-7 
ur =1600 

lc=160/ 100 
lg=0 . 8/1000 
A=5/ 10000 
N=1200 

Rc=lc/ (uo*ur*A) 
Rg  =  lg / ( uo  *  A ) 
R=Rc+Rg 

L=N*N/R 
disp (L) 


Scilab  code  Exa  7.15  Probleml5 

N  =  1 00 
11=0.15 
12=0 . 3 
13=0.45 
A  =  0 . 001 

url =1447 
ur2  =5969 
ur3  =47750 

f lux=0 . 6E-3 

B=f lux/A 

HI =B/ ( uo  *ur 1 ) 

H2  =  B  /  (  uo  *  ur  2  ) 
H3=B/(uo*ur3) 
disp (HI , H2 , H3) 
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F=H1*11+H2*12+H3*13 
disp (F) 

I  =  F/N 
di sp ( I ) 


Scilab  code  Exa  7.16  Problcml6 

PI 1 =1500 
f  1=50 
P12=3000 
f  2=75 


Y=  [1  ,  f  1  ;  1  ,  f  2] 

X= inv (Y) * [Pll/fl ; P12/f2] 

Phl=X (1) *f  1 
Pe 1 =X  (2) *f l*f 1 
Ph2  =  X  (1) *f 2 
Pe2  =  X  (2) *f 2*f 2 

disp (Phi , Pel , Ph2 , Pe2) 


Scilab  code  Exa  7.17  Problcml7 


uo  =  (4*°/„pi)  +1E-7 
F=750* 1 
disp (F) 

dm =(16+20) /200 
lc  =  dm*°/0pi 
Hc=F/lc 
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disp (He ) 


f lux=l . 25/1000 
Ac=(20-16) *2 . 5/10000 
Bc=f lux/ Ac 
disp (Be) 

R=F/ flux 
disp (R) 

u=Bc/Hc 
di sp (u) 

ur =u/uo 

disp(ur)  //calculation  mistake  in  the  book 

here 


Scilab  code  Exa  7.18  Probleml8 

f luxe  =  1 . 5E-3 
Ac=12/10000 
Bc  =  f luxe /Ac 
He  =  1 250 
lc  =  0 . 3 
F  =  Hc  *  1 c 


Fa=200* 1 
Fb=600*0 . 75 

Fc=F-Fb+Fa 
Nc=Fc / 0 . 5 

disp(Nc)  // // clockwise5 
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Scilab  code  Exa  7.19  Problcml9 


uo  =  (4*y„pi)  *lE-7 

ur  =  2000 

V  =  200 

f  =50 

N=1600 

Ac  =  5/ 10000 

f lux_max=V/ (4 . 44* f *N) 
Braax=flux_max/Ac 
di sp ( Bmax ) 

Rc=20/100/(uo*ur*Ac) 

Rg=0 . 5/1000/ ( uo  *  Ac ) 
R=Rc+Rg 

imax=flux_max*R/N 
di  sp  (  iraax ) 

Wfmax=l/2*R*(flux_max~2) 
disp (Wf max) 

percent =Rg/R* 100 
disp (percent ) 
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Chapter  8 
Transformers 


Scilab  code  Exa  8.1  Probleml 

VI =3300 
f  =50 
N1 = 100 
N2  =  300 
Z=100+35*%i 

f  lux_max  =  Vl/(sqrt  (2)  *°/„pi  *f  *N1 ) 

V2=V1*N2/N1 

I2=V2/Z 

I1=N2/N1*I2 

S  =  V1 *  con j  (II) 

Z1=Z* (N1/N2) “ 2 

di sp ( f lux_max ) 
disp(Il  ,  12) 
disp(real(S) , imag ( S ) ) 
disp(Zl) 


Scilab  code  Exa  8.2  Problem2 
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a=2400/240 


R=0.2+a*a*2/1000 
X=0 . 6+a*a*6/ 1000 
disp (R , X) 

Rlv=l/ a/ a*R 
Xlv=l/ a/ a*X 
disp (Rlv , Xlv) 

12=150*1000/2400 
Z  =  R  +  %  i  *  X 
V= 12  *Z 

pcnt=norm(V) /2400*100 
di sp (V , pent ) 

Im=2400/ 1600 
Ii =2400/ 10000 
Io  =  Ii-°/„i*Im 

pf=cos (atan(imag(Io)/ real (Io) ) ) 
disp(Io ,pf) 


Scilab  code  Exa  8.3  Problcm3 

P=200E3 
Vl=11000 
V2  =4 1 5 
f  =50 
N2  =80 

N1=V1/V2*N2 
a=Ml/N2 
di sp ( a) 

I2=P/V2 
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disp(I2) 

I 1 = 12 / a 
disp(Il) 

Z2  =  V2  * V2/P 
disp(Z2) 
Z21=Z2*a*a 
disp(Z21) 


Scilab  code  Exa  8.4  Problem4 

VI =200 
f  =50 
N1 = 150 

Ac  =  1 0  *  5 / 10000 


Flux_max  =  Vl/(sqrt  (2)  *°/„pi  *f  *N1 ) 

Braax=Flux_max/Ac 

Hraax  =250 

1=2* (25+10) +2* (20+10) 

1=1/100 
AT_max=Hmax*l 
Im_max=AT_max/150 
Im_rms = Im_max/ sqrt (2) 
disp(Im_rms) 

Fe.loss =23000 

Cv  =  2* (25+2*10) *10*5  +  2*20*10*5 
Cv=Cv/ 1000000 
Cl=Fe_loss*Cv 
di sp ( Cl ) 

Ii=150/Vl 


129 


25 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

1 

2 

3 

4 

5 

6 

7 

8 


disp(Ii) 


Scilab  code  Exa  8.5  Problcm5 

function  [x,y]=polar_to_cart (r, theta) 
theta  =  theta/180*°/0pi 
x  =  r *  cos (theta ) 
y=r*sin (theta) 
endf unction 

[Zr,Zc]=polar_to_cart (5,30) 
Z2=complex(Zr,Zc) 

V  =  200 
N2  =75 
N1 =100 

I2=V/Z2 

I1=N2/N1*I2 

disp(norm(Il)) 

disp(”  lagging”  ,  real ( I 1 ) /norm ( I 1 ) ) 


Scilab  code  Exa  8.6  Problcm6 


R  =  60 
V  =  6 

Rs  =2400 

a=  sqrt ( Rs /R) 
di sp ( a) 

Pl_max=l/2*(R*R)/ (Rs+Rs) 
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Il=V/2/Rs 


11  =  1 1  *  a 
disp ( II ) 
Vl=V/2/a 
disp (VI ) 


Scilab  code  Exa  8.7  Problcm7 


f  =50 

a=1100/220 
Rhv  =  0 . 125 
Xhv  =  0 . 625 
R1 v  =  0 . 005 
XI v  =  0 . 025 

Zhv  =  Rhv  +  °/0i  *  Xhv 
Zlv  =  Rlv  +  °/0i*Xlv 

Zl=Zhv+a*a*Zlv 
di sp ( Z 1 ) 

Z2=Z1 /a/a 
di sp ( Z2 ) 

Zpul =Z1 *0 . 05/1 .1/1.1 
di sp ( Zpul ) 

Zpu2=Z2*0 .05/0.22/0.22 
di sp ( Zpu2  ) 


Scilab  code  Exa  8.8  Problem8 

P=600E3 

a=2400/600 
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r 1 =0 . 05 
r2  =0 . 004 
x 1 =0 . 025 
x2  =0 .016 
R1 =1667 
Xm=417 

/////As  seen  from  the  LV  side 

Zlv  =  r2  +  %i*x2  +  l/a/a*(rl  +  °/0i*xl) 
disp(Zlv) 

RiLV=Rl /a/a 
disp(RiLV) 

XmLV=Xm/ a/a 
di sp ( XmLV ) 


Zpu=Zlv*0.6/0.6/0.6 
di sp ( Zpu ) 

Ri  =  RiLV*0 . 6/0 . 6/0 . 6 
disp (Ri ) 

Xm=XmLV*0 . 6/0 . 6/0 . 6 
disp (Xm) 


Scilab  code  Exa  8.9  Problcm9 

P  =  50  E3 
a=2200/220 

II  III OC  Parameters 

Poc  =405 
Ioc  =5 
Voc  =220 
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/////SC  Parameters 

Psc  =805 
Isc=20 . 2 
Vsc  =95 


Y0=loc/Voc 
Gi=Poc/Voc/Voc 
Bm=sqrt (Y0*Y0-Gi*Gi) 

Z=Vsc/Isc 
R=Psc/ Isc/ Isc 
X  =  sqrt (Z*Z-R*R) 

/////////// Referred  to  HV  side 

GiHV=Gi/a/a 

disp(GiHV) 

BraHV=Bm/ a/ a 
di sp ( BraHV ) 
disp (R) 
disp (X) 

//////////// Referred  to  LV  side 

disp(Gi) 
disp  (Bin) 

RLV=R/a/a 
di sp ( RLV ) 

XLV=X/ a/a 
disp (XLV) 


////////// per  unit 

GiPU=GiHV/0 .0103 
BmPU=BmHV/0 .0103 
RPU=R/96 . 8 
XPU=X / 96 . 8 
disp (GiPU) 
di sp ( BmPU ) 
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disp (RPU) 
disp (XPU) 


Scilab  code  Exa  8.10  ProblemlO 

P  =  50  E3 
a=2200/110 

/////OC  Parameters 

Poc  =400 
Ioc  =  10 
Voc  =  1 10 

Y0=loc/Voc 

Gi=Poc/Voc/Voc 

disp(Gi) 

Bm=sqrt (Y0*Y0-Gi*Gi) 
disp  (Bin) 

/////////// Referred  to  HV  side 

GiHV=Gi/a/a 

disp(GiHV) 

BraHV=Bm/ a/a 
di sp ( BraHV ) 

////////// per  unit 

GiPU=Gi*P/Voc/Voc 
BmPU  =  Bm  *  P  /  Vo  c  /  Vo  c 
disp (GiPU) 
di sp ( BraPU ) 


Scilab  code  Exa  8.11  Problcmll 
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P=25000 
VI =2200 
V2  =  220 
a=Vl/V2 
f  =50 
r  1  =2 

r2  =0 . 025 
x  1  =7 
x2  =0 . 07 
Xm  =  1 6000 

Im=V2/Xm 

disp(”pf  =  0,  90  degree  lag  Po=0”  ,  Im) 


R=rl+a*a*r2 
X=xl+a*a*x2 
Z  =  sqrt (R*R+X*X) 

If 1=P/V1 
Vsc  =  Z* If  1 

Vsc_percent=Vsc/Vl*100 
disp(Vsc ,Vsc_percent) 
Im=Vsc/Xra 

Im_percent=Im/Ifl*100 
dispClm, Im_percent) 

Sc_pf=cos (atan(X/R) ) 
disp (Sc_pf ) 

// ////////// per  unit 

ZBpu  =  10  00*2 . 2*2 . 2/25 
r lpu=rl /193 . 6 
r2pu=a*a*r2/193.6 
xlpu=xl / 193 . 6 
x2pu=a*a*x2/193 . 6 

Rpu=rlpu+r2pu 

Xpu=xlpu+x2pu 
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Xrapu=Xm/ 193.6 
disp (Rpu , Xpu , Xmpu) 


Scilab  code  Exa  8.12  Probleml2 

P=10000 
vl =2300 
v2  =  230 
a=vl /v2 
f  =50 
r 1 =3 . 96 
r2  =  0 . 0396 
xl =15 . 8 
x2  =0 . 158 

R=rl+a*a*r2 

X=xl+a*a*x2 

V2=v2*a 

I=P/V2 

theta=acos (0.80) 

Vl=V2+I*(R*cos (theta)+X*sin(theta) ) 
disp (VI ) 

VR= (Vl-V2)/V2 
disp (VR) 

pf =X/sqrt (R*R+X*X) 
theta2=acos (pf) 

Il=I*(cos  (theta)  - °/„i*sin(theta)  ) 
Ic=real (Il*tan(theta2) ) -imag(Il) 
Rating_Cap=V2*Ic 
disp(Rating_Cap) 

VI  =  V2 
disp (VI ) 
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Scilab  code  Exa  8.13  Problcml3 


P=100000 
N 1 =400 
N2  =  100 
a=Nl/N2 
r 1 =0 . 3 
r2  =0 .015 
xl=l . 1 
x2  =0 . 055 
VI =2400 

R=rl+a*a*r2 
disp (R) 

X=xl+a*a*x2 
disp (X) 

I1=P/V1 
pf =0 . 8 

theta=  acos (pf ) 

Vd=Il*(R*cos (theta)+X*sin(theta) ) 

V R=Vd/Vl *100 
V2= ( VI -Vd) / a 
disp (VR , V2 ) 

Vd=Il*(R*cos (theta) -X*sin(theta) ) 

V R=Vd/Vl *100 
V2= ( VI -Vd) / a 
disp  (VR  ,  V2 ) 

pf = cos ( atan (R/X) )  /////wrong  in  the  book 

disp  (”  leading”  ,  pf  ) 
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Scilab  code  Exa  8.14  Probleml4 


P=10000 
vl =2300 
v2  =  230 
a=vl / v2 
f  =50 
r 1 =3 . 96 
r2  =  0 . 0396 
xl  =  15 . 8 
x2  =0 . 158 

R=rl+a*a*r2 
X=xl+a*a*x2 
V2  =  v2  *  a 
I=P/V2 
pf =0 . 8 

theta  =  acos  (pf ) 

Vl=V2+I*(R*cos (theta)+X*sin(theta) ) 

Pi=75*Vl*Vl/V2/V2 

Pc  =  I  *  I *R 

Pl=Pi+Pc 

P0=P*pf 

effi=P0/ (P0+P1) *100 
disp ( ef f i  ) 

VI  =  V2 
Pi  =75 

pf 2=X/sqrt (R*R+X*X) 

I=P0/Vl/pf 2 
Pc  =  I  *  I *R 
Pl=Pi+Pc 
P0=P*pf 

effi=P0/ (P0+P1) *100 
disp ( ef f i  ) 

I  =  sqrt ( Pi /R) 

Load=V2*I 


138 


37 

38 

39 

40 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 


P0=Load*pf 
PI =2*  Pi 

effi_max=P0/ (P0+P1) *100 
di sp ( ef  f i_max ) 


Scilab  code  Exa  8.15  Probleml5 

P  =  1 5  E3 
a=3000/250 

/////OC  Parameters 

Poc  =  105 
Ioc  =0 . 62 
Voc  =250 


/////SC  Parameters 

Psc  =360 
Isc=5 . 2 
Vsc  =  157 


Y0=loc/Voc 
Gi=Poc/Voc/Voc 
Bm=sqrt (Y0*Y0-Gi*Gi) 

Z=Vsc/ Isc 
R=Psc/ Isc/ Isc 
X  =  sqrt (Z*Z-R*R) 

///////////  Referred  to  HV  side 

GiHV=Gi/a/a 

disp(GiHV) 

BraHV=Bm/ a/a 
di sp ( BmHV ) 
disp (R) 
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disp (X) 


////////// per  unit 

ZB  =  3 . 12/5 . 2 
GiPU=Gi*ZB 
BmPU=Y0*ZB 
RPU=R/ZB/1000 
XPU=X/ZB/1000 
disp (GiPU) 
di sp ( BmPU ) 
disp (RPU) 
disp (XPU) 


pf =0 . 8 

theta=  acos (pf ) 

V2=250*a 

I=P/V2 

Vd=I*(R*cos (theta) -X*sin(theta) ) 

VR=-Vd/V2*100 

disp (VR) 

Pi = 105 
Pc  =  I  *  I *R 
Pl=Pi+Pc 
P0=P*pf 

effi=P0/ (P1+P0) *100 

pf2=cos (atan(R/X) ) 
disp  (”  leading”  ,  pf2) 

IPU= sqrt (Pi/P/RPU) 
ef f i_max= IPU 
di sp ( ef  f i_max ) 
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Scilab  code  Exa  8.16  Problcml6 


vl =2500 
v2  =  250 
P=25000 
Pc  =  130 
Pcf 1=320 
pf =0 . 8 

V2  =2750 

P0=P*pf 
Pl=Pcf 1+Pc 
effi=P0/ (P0+P1) *100 
disp ( ef f i  ) 

I2=P/v2 
I 1 = I2+P/ vl 
kVA=Il*vl/1000 
P0=kVA*pf *1000 
effi=P0/ (P0+P1) *100 
disp ( ef f i  ) 


Scilab  code  Exa  8.17  Probleml7 


a  =  10 

Load=120000 
V  =  400 

I=Load/V/ sqrt (3) 
Vll=a*V/ sqrt (3) 
111= sqrt (3) *I/a 
all=Vll/V 
disp (Vll , Ill , all) 
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Vll=a*V*sqrt (3) 
Ill=I/a/ sqrt (3) 
all=Vll/V 
disp (Vll , Ill , all) 


Scilab  code  Exa  8.18  Probleml8 


f  =50 
N 1 =500 


Pin=60 
Io  =  0 . 4 
Vin=220 
r  =  0 . 8 


Pci=Io*Io*r 
Pi=Pin-Pci 
disp (Pi ) 

theta=acos (Pin/Vin/Io) 
Im=Io*sin(theta) 
Xm=Vin/Ira 
disp (Xm) 

Iio=Io*cos (theta) 

Ri=Vin/Iio 

disp (Ri ) 


Scilab  code  Exa  8.19  Problcml9 

P  =  1 5  E3 
vl =2200 
v2  =  220 
a=vl / v2 
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/////OC  Parameters 

Poc  =  185 
Ioc  =2 . 72 
Voc  =220 


/////SC  Parameters 

Psc=197 
Isc=6 . 3 
Vs  c  =  1 1 2 


Pi=Poc 

disp(Pi) 

IHVf l=P/vl 

Pcf l=IHVf l*IHVf l/Isc/Isc*Psc 
disp(Pcfl) 

pf =0.85 
Po=P*pf 
Pl=Pi+Pcf 1 
effi=Po/ (Po+Pl) *100 
disp ( ef f i  ) 

Z=Vsc/ Isc 
R=Psc/ Isc/  Isc 
X  =  sqrt (Z*Z-R*R) 
pf =0 . 8 

theta  =  acos  (pf ) 

Vdl=IHVfl*(R*cos (theta)+X*sin(theta) ) 
Vd2=IHVfl*(R*cos (theta) -X*sin(theta) ) 
VRl=Vdl/vl *100 
VR2=Vd2/vl *100 
disp  ( VR1  ,  VR2  ) 


Scilab  code  Exa  8.20  Problem20 
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P=50000 
lr =0 . 9 

P0=P*1*0 . 9 

eff i=0 . 974 

Pl=(l-effi)/effi*PO 

Pi=Pl/2 

Pcf 1 =Pi / lr / lr 

pf =0 . 8 
P0=P*pf 
Pl=Pi+Pcf 1 
effi=P0/ (P0+P1) *100 
disp ( ef f i  ) 

P0=P/2* lr 
Pl=Pi+Pcf 1/2/2 
eff i=PO/(PO+Pl) *100 
disp ( eff i  ) 

/////////calculation  mistakes  in  the  book 


Scilab  code  Exa  8.21  Problem21 

P=500E3 
ef f i=0 . 95 

A  = [ 1  ,  1 ; 1  ,0.6*0. 6] 

Pa= inv (A)*[P*(l-effi)/effi;P*0.6*(l-effi)/effi] 
Pi  =  Pa  (1) 
disp (Pi ) 

Pc=Pa (2) 
disp(Pc) 

Pl  =  Pi +0 . 75*0 . 75*  Pc 
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effi=P*0 . 75/ (P*0 . 75+P1) 
disp ( ef f i  ) 


Scilab  code  Exa  8.22  Problem22 

vl =2200 
v2  =  220 
f  =50 

emf perturn=12 

Nl=floor(vl/emfperturn) 
N2=floor (v2/ emfperturn) 
disp(Nl) 
disp (N2) 

fluxmax=emfperturn/4.44/f 
Braax=l . 5 
Ac=f luxmax/Bmax 
disp(Ac) 


Scilab  code  Exa  8.23  Problem23 

vl =3300 
v2  =  600 
f  =50 

Ac=25/10000 

1  =  1.2 


Braax=l . 2 

fluxmaxl=Bmax*Ac 
Nl=vl/4.44/f/fluxmaxl 
N2=v2/4 . 44/ f / f luxmax 1 
disp ( N 1 , N2 ) 
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12=20 
pf =0 . 8 
I1=N2/N1*I2 

Il  =  Il*(pf-0/0i*sin(acos  (0.8)  )  ) 
di sp ( 1 1 ) 

Hraax  =250 
ATmax=Hmax*l 
immax=ATmax/Nl/sqrt (2) 

Ii  =0 

1 1  =  - °/„ i  *  immax  +  Il 

disp  (”  lagging”  ,  real ( 1 1 ) /norm ( 1 1 )  ,  norm(Il)) 


Scilab  code  Exa  8.24  Problem24 

P=50000 
vl =2400 
v2  =  240 
f  =50 
a=vl / v2 
Pd=375 
pf =0.4 

I0=Pd/vl/pf 
I0_prime=10*a 
disp ( I0_prime  ,  pf ) 


Scilab  code  Exa  8.25  Problem25 

vl =220 
v2  =  1 1 0 

zl=0.32  +  °/0i*0.85 
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z2  =0 . 1 1  +  % i *  0 • 27 
a=vl / v2 

z=zl+a*a*z2 
IHV=vl/ norm (z) 
ILV= IHV*a 
disp (IHV , ILV ) 


Scilab  code  Exa  8.26  Problem26 

P=1000000 
vl =1 1000 
v2  =  230 
f  =50 
Vsc=310 
Psc=5210 

Isc=P/ vl 
Z=Vsc/ Isc 
R=Psc/ Isc/ Isc 
X  =  sqrt (Z*Z-R*R) 


pf  =  1 

theta=  acos (pf ) 

Vd=Isc*(R*cos (theta)+X*sin(theta) ) 
VR  =  Vd/ vl *  100 
Vl=vl+Vd 
disp (VR , Vl) 

pf =0 . 8 

theta  =  acos  (pf ) 

Vd=Isc*(R*cos (theta)+X*sin(theta) ) 
VR  =  Vd/ vl *  100 
Vl=vl+Vd 
disp (VR , Vl) 
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pf =0 . 8 

theta=  acos (pf ) 

Vd=Isc*(R*cos (theta) -X*sin(theta) ) 
mistake  in  the  book  at  this  point 
VR  =  Vd/ vl *  100 
Vl=vl+Vd 
disp (VR , VI) 


Scilab  code  Exa  8.27  Problem27 

P=10000 

vl =2200 

v2  =  220 

r  1  =4 

r2  =0 . 04 

x  1  =5 

x2  =0 . 05 

a=vl /v2 

R=rl+a*a*r2 

X=xl+a*a*x2 

Il=P/vl 

pf =0 . 8 

V2  =  v2  *  a 

thet a= acos (pf  ) 

Vd=Il*(R*cos (theta)+X*sin(theta)  ) 
VR  =  Vd/ vl *  100 
Vl=vl+Vd 
disp (VR , VI) 

pf=cos ( at an (X/R) ) 
disp  (”  leading”  ,  pf  ) 


Pr=P*pf 

Q=-P*sin(acos (pf) ) 
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disp (Q , Pr ) 


Scilab  code  Exa  8.28  Problem28 

P=20000 
vl =200 
v2  =400 
f  =50 
VI =600 
V2  =  200 

a=Vl/V2 
di sp ( a) 

Il=P/v2 
VA=V1 *11 
disp (VA) 

di sp ( 1 1 ) 

12 =30000/ V2 
I s  =  12  - 1 1 
di sp ( 12 ) 
disp(Is) 

VAtr ans  =  v2  *  1 1 
VAcond =30000 -VAtr an s 
disp(VAtrans) 
disp (VAcond) 


Scilab  code  Exa  8.29  Problcm29 

P=100000 
vl =1 1000 
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v2  =400 


VI =400 

Vp=Vl/sqrt (3) 

a=ceil (vl/Vp) 
di sp ( a) 

I1  =  P/ sqrt (3)  /VI 
Ip  =  Il 

disp (II , Ip) 

Il=P/sqrt (3)/vl 
Ip=Il/sqrt (3) 
disp ( II , Ip) 


150 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 


Chapter  9 

EMF  and  Torque  in  Electric 
Machines 


Scilab  code  Exa  9.1  Problcml 


f  =50 

f lux=0 .016 
S  =36 
P  =  6 
N  =  1 0 


m  =  S/P 

gammaa  =  °/„pi  /  m 

Kb=sin (m* gammaa/2) /m/ sin (gammaa/2) 

Nph  =  S  *  N  *  2/2 

Ep=4.44*Kb*f*Nph*flux 

kVA=Ep*N/1000 

disp(kVA) 

ra=S/2/P 

gammaa  =  7„pi  /  2/m 

Kb=sin (m* gammaa/2) /m/ sin (gammaa/2) 
Nph=S*N*2/2/2 
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Ep=4.44*Kb*f*Nph*flux 
El ine =Ep * sqrt (2) 
kVA=Ep*N*2/1000 
disp(kVA) 

m=S/3/P 

gammaa  =  y„pi  /  3/m 

Kb  =  sin (m* gammaa/2) /m/ sin ( gammaa/2) 

Nph=S*N*2/2/3 

Ep=4.44*Kb*f*Nph*flux 

El ine =Ep * sqrt (3) 

kVA=Ep*N*3/1000 

di sp ( kVA ) 


Scilab  code  Exa  9.2  Problem2 

S  =54 
P  =  6 

m=S/3/P 

gammaa  =  °/„pi  /3/m 

Kbl  =  s in (m*  gamma a /2) /m/ sin ( gammaa/2) 
Kb3=sin (m* gammaa/2 *3) /m/ sin (gammaa/2*3) 
Kb5=sin (m*gammaa/2*5) /m/ sin (gammaa/2*5) 

disp (Kbl ) 
di sp ( Kb3 ) 
di sp ( Kb5 ) 


Scilab  code  Exa  9.3  Problcm3 


f  =50 
n  =  500 
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ra  =  5 
N  =  1 2 

f lux=0 . 025 

P=120*f/n 
S=m*3*P 
Nph=S*N*2/2/3 
gamraaa  =  7„pi  /  3/m 

Kb  =  sin (m* gammaa/2) /m/ sin ( gammaa/2) 

polepitch=S/N 

coilpitch=13 

spa=(polepitch-coilpitch)*  gamma a 
Kp  =  cos ( spa/2) 

Ep=4.44*Kb*Kp*f*Nph*flux 

di sp ( Ep ) 

Eline=sqrt (3) *Ep 
disp(Eline) 


Scilab  code  Exa  9.4  Problem4 


f  =50 
P  =  6 

ns=120*f /P 
disp (ns ) 
wm=2*%pi/60*ns 
disp (wm) 


Scilab  code  Exa  9.5  Problem5 

F2  =  850 
F 1 =400 

lambda  =  123 .6/180* %pi 
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ppp  =  l  .408/1000 
theta=%pi - lambda 

Fr=sqrt (Fl*Fl+F2*F2-2*Fl*F2*cos (theta) ) 

f  luxr  =  2/°/„pi*ppp*Fr 
di sp ( f luxr ) 


Scilab  code  Exa  9.6  Problcm6 


P  =  6 

N=1000 

f =P*N/120 
ns  =  120*  f /4 
s  =  0 . 04 
n=(l-s) *ns 
disp (n) 


Scilab  code  Exa  9.7  Problcm7 


P  =  4 
f  =50 


ns=120*f /P 
nl = - 1500 
sl=(ns-nl)/ns 
f 2=sl*f 
SRV  =  s 1*80 
disp (f 2 , SRV ) 


154 


12 

13 

14 

15 

16 

17 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

1 

2 

3 

4 


n2=1000 
s2  =  (ns -n2 ) /ns 
f 2=s2*f 
SRV=s2*80 
disp (f 2 , SRV ) 


Scilab  code  Exa  9.8  Problcm8 


betaa=acos (0) *2/5 
coilpitch=%pi-betaa 
disp (coilpitch/%pi *180) 


P  =  6 
S  =72 
ra  =  S/P 

gammaa  =  7„pi  /  m 
cpis  =  150/180*°/„pi 
betaa=%pi-cpis 
Kpl=cos (betaa/2) 
Kp3=cos (5* betaa/2) 
Kpl3=cos (13*betaa/2) 
disp (Kpl ) 
di sp ( Kp3 ) 
disp(Kpl3) 


Scilab  code  Exa  9.9  Problcm9 


P  =  2 
f  =50 
S  =42 
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ra=S/3/P 

gammaa  =  °/„pi  /3/m 

Kb  =  sin (m* gammaa/2) /m/ sin ( gammaa/2) 

coilpitch=17 

polepitch=S/2 

spa=(polepitch-coilpitch)*  gamma a 
Kp  =  cos ( spa/2) 


N  =  S  *2 
Nc  =  N/2 

Coilsperphase=Nc/3 

Nph=Coilsperphase/2 

Eline=2300 

flux=Eline/sqrt (3) /4 . 44/Kb/ Kp/f/Nph 
di sp ( flux ) 


Scilab  code  Exa  9.10  ProblcmlO 


f  =50 
ns=1000 
P=120*f /ns 
disp (P) 
s= (ns -940) / ns 
disp ( s  *  100 ) 
nrr =ns -940 
disp (nrr) 
nrs=1000 
disp (nrs ) 
s=2*  s 

n=1000* (1 -s) 
disp (n) 


Scilab  code  Exa  9.11  Problemll 
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ns=1000 
s= (ns -940) / ns 
news  =  2* s 
n=1000*(l-news) 
disp (n) 


Scilab  code  Exa  9.12  Problcml2 


P  =  8 
f  1=60 
f  2  =50 
If  1=5 

If 2=f l/f2*If 1 
disp(If2) 
nA2  =  120*f2/P 
disp (nA2 ) 


Scilab  code  Exa  9.13  Probleml3 

f  1=60 
f  2  =50 
PI  =  12 

P2=f 2/f 1*P1 
disp (P2) 
nset =120*f 1/PI 
disp(nset) 


Scilab  code  Exa  9.14  Probleml4 


ra  =  6 
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gammaa  =  y„pi  /  m 

Kb = sin (m*gammaa/2) /m/ sin (gammaa/2) 


m  =  4 

Kb  =  sin (m*  gammaa/2) /m/ sin ( gammaa/2) 


Scilab  code  Exa  9.15  Problcml5 


f  =50 
n  =  965 

P=f loor (120*f/n) 
disp (P) 

s= ( 1000 -965) /1000 
di sp ( s ) 
f 2=s*f 
disp (f 2) 

nsr=1000-965 
disp (nsr ) 
nrr  =  120*  f 2/P 
nsrf=1000-965-nrr 
di sp ( nsrf ) 


Scilab  code  Exa  9.16  ProblcmlG 


f  =50 
P  =  4 
a  =  2 


ns=120*f /P 
n= 1440 
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7  s=l-n/ ns 

8  di sp ( s ) 

9  f 2=s*f 

10  disp  (f  2) 

11  Elp=400 

12  E2f  =400/  a 

13  E2f  2  =  E2f  *f  2/f 

14  E21=sqrt (3) *E2f 2 

15  disp(E21) 
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Chapter  10 
DC  Machines 


Scilab  code  Exa  10.1  Problcml 


P  =  6 
S  =36 

coilside=2 
N  =  8 
A  =  2 

d=0 . 25 
1=0 . 18 
Ia  =  10 

avgf lux  =  0 . 8 
n=1200 


Z  =  S*  coilside  *N 

f lux  =  %pi *d/P*l*  avgf lux 

Ea=flux*n*Z/60*P/A 

Pm=Ea* la 

disp (Pm) 

T  =  Pm/2  /  °/„pi  /  n*60 
disp (T) 
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Scilab  code  Exa  10.2  Problem2 


P  =  4 
A  =  4 

ra=0 . 145 
1=0.21 
Z=2*33* 1 1 
K  =  Z*P/2/°/„pi/A 
disp (K) 


Ap  =  2*  %pi *  ra/P  *0 . 7* 1 
Bar c  =0 . 8 
f lux=Ap*Barc 
n=1200 

Ea  =  K*flux*2*°/0pi*n/60 
disp (Ea) 

Ia=240 
Ic= Ia/A 
disp(Ic) 

T=K*f lux*Ia 
disp (T) 

Pg=Ea* la 
disp (Pg) 


Scilab  code  Exa  10.3  Problcm3 

P=5000 
Vt  =215 
n=1000 
Ra=0 . 4 
Ia=P/Vt 
Eag=Vt+Ra* la 
Eara=Vt -Ra* la 
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newn=Eam/Eag*n/l . 1 
di sp ( newn ) 


Scilab  code  Exa  10.4  Problem4 

P=200000 
V  =  400 
n  =  600 
Z  =  864 
PI =8000 

Ia=P/V 
Ra=Pl/Ia/Ia 
Ea=V- Ia*Ra 
flux=Ea*60/n/Z 
di sp ( flux ) 


Scilab  code  Exa  10.5  Problem5 
N=1800 

///////  from  figure 

Voc  =250 
If =4 . 1 
Rf =Voc/If 
disp(If) 
disp (Rf ) 

Rfcrit=150/2 

disp(Rfcrit) 

V3=120 

Ncrit=N*V3/ 150 
disp (Ncrit ) 
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Rf 220  =  220/3 . 2 
Rext  =  Rf 220 -Rf 
di sp ( Rext ) 


Scilab  code  Exa  10.6  Problcm6 


P  =  4 
V  =  230 
Z  =  888 
Ra=0 . 8 
f lux=5 . 4E-3 
A  =  2 


11=2 
If =0 . 6 
Ia= II - If 
Ea=V- Ia*Ra 
n0=Ea*60/ flux/Z/P*A 
disp (nO) 

T  =  29 . 6 

Ia  =  T*2*°/„pi/flux/Z/P*A 

Il=Ia+If 

disp ( II ) 

Ea=V-Ra* la 
n=Ea*60/flux/Z/P*A 
disp (n) 

sr  =  ( 1 -n/nO ) *  100 
disp(sr) 


Scilab  code  Exa  10.7  Problem7 
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nl =1000 
V  =  230 
Ia  =  75 
Ra=0 . 1 
Eal =V- Ia*Ra 
Rf =275 
If =V/Rf 
Ke=Eal/If /nl 

n2=1200 
Ia= 125 
Ea2=V-Ia*Ra 
If 2=Ea2/Ke/n2 
Rf 2=V/If 2 
Rf  ext  =  Rf 2 -Rf 
disp(Rfext) 


Scilab  code  Exa  10.8  Problem8 

V  =  1 15 
Ial =25 
Ra=0 . 3 

nl = 1450 
Eal =V- I al *Ra 
Ke=Eal /nl 

n2=1200 

Ea2  =  Ke  *n2 

Ia2=3/ 4* Ial 

Raext=(V-Ea2)/Ia2-Ra 

di sp ( Raext ) 

ef  f  ia  =  Ea2/V*  100  ///calculation  mistake  in  the 

book  at  this  point 

disp ( ef f ia) 
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V=Ea2+Ia2*Ra 
effia=Ea2/V*100 
disp ( ef f ia) 


Scilab  code  Exa  10.9  Problem9 

/////solving  quadratic 

Ra=0 . 5 
P=8000 
V  =  230 
Ea  =  V 
n=1200 
Ke=V/n 

p  =  poly ( [P*Ra , -V  ,  1]  ,  ”w”  coeff”) 

w=roots (p) 

Ea  =  w(l)  ///rejecting  small  value 

n=Ea/Ke 
disp (n) 

T  =  P/  2/yopi/n*60 
disp (T) 

Ia= (V-Ea) /Ra 
disp ( Ia) 

Kt=T/Ia 


p  =  poly ( [-Kt*V/Ra , Kt*Ke/Ra  ,0 . 6E-4]  ,  ”w”  coeff  ”  ) 
w=roots (p) 
n=w (2) 
disp (n) 


Scilab  code  Exa  10.10  ProblcmlO 
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V  =  300 
Ea  =  V 
nl =1200 
n2  =  1 100 

Kaphi  =  Ea/2/°/„pi/nl*60 
T  =  350 

Ia=T/Kaphi 
disp ( la) 

Ea=V*n2/ nl 
Pm=Ea* la 
disp (Pm) 

Ra=(V-Ea) /la 
disp (Ra) 

disp (T) 

Ea=600 - Ia*Ra 

n  =  Ea*60/Kaphi/2/°/0pi 

Pm=Ea* la 

disp (Pm) 

di sp (n) 


Scilab  code  Exa  10.11  Problemll 

V  =  300 
Ea  =  V 
nl =1200 
n2  =  1 100 

Kaphi  =  Ea/2/°/„pi/nl*60 
T  =  350 

Ia=T/Kaphi 
Ea=V*n2/ nl 
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Ra=(V-Ea) /la 


Kaphi =Kaphi / 2 

n  =  Ea/Kaphi*60/2/°/0pi 

Pm=Ea* la 

T  =  Kaphi *  la 

disp (T) 

disp (Pm) 

di sp (n) 


Scilab  code  Exa  10.12  Probleml2 


V  =  600 
Ia=40 
R  =  0 . 5 
Ea=¥- Ia*R 


n  =  500 

Ka  =  Ea/Ia/2/°/„pi/n*60 
T=Ka* la* la 


n2  =450 

T  =  T*n2  *n2 /n/n 
Ia=sqrt (T/Ka) 

Ea  =  Ka*Ia*2*°/„pi*n2/60 
Rtotal=(¥-Ea)/Ia 
Rext  =  Rtotal -R 
di sp ( Rext ) 


Scilab  code  Exa  10.13  Probleml3 

¥  =  220 
Ra  =  l 
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Rse  =0 . 4 
Ial =20 

Ia2=sqrt(Ial*Ial*0.7*0.7*0.7) 

Eal=V-Ial*(Ra+Rse) 

Ea2  =  Ia2  *0 . 7*Eal/Ial 
Rext=(V-Ea2)/Ia2-Ra-Rse 
di sp ( Rext ) 


Scilab  code  Exa  10.14  Probleml4 

V  =  250 
Ia  =  25 
n=1000 

KaNfNse  =  V/Ia*60/2/°/0pi/n 
T=KaNfNse*Ia*Ia 

Ia=sqrt  (V*2/KaNfNse/  sqrt  (KaNfNse*2*°/0pi*n/60*2*%pi*n 
/60/T/2) ) 

w=sqrt  (KaNfNse*2*°/0pi*n/60*2*°/0pi*n/60/T/2)  *  la 
n=w*60/2/%pi 

disp (n) 
disp ( la) 


Scilab  code  Exa  10.15  Probleml5 

V  =  230 
Rf =120 
Ra=0 . 15 
If =V/Rf 
Psh=V* V/Rf 

1=14 . 5 
Pin  =  I  * V 
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Ia=I - If 

Pk=Pin-Ia*Ia*Ra 

1=215 
Ia=I - If 

Pl= la* Ia*Ra+Pk 
Pin=V* I 

effi=(Pin-Pl)/Pin 
disp ( ef f i  ) 

Ia=sqrt (Pk/Ra) 
Il=Ia+If 
PI =2*  Pk 
Pin=V* II 

effi=(Pin-Pl)/Pin 
disp ( ef f i  ) 


Scilab  code  Exa  10.16  Probleml6 

V  =  250 
Rf =125 
Ra=0 . 2 
If =V/Rf 
1  =  16 

IaO=I - If 

Pk=¥*IaO-IaO*IaO*Ra+¥*If 

1  =  152 
Ia=I - If 

Pl= la* Ia*Ra+Pk 
Pin=¥* I 

effi=(Pin-Pl)/Pin 
disp ( ef f i  ) 

11=152 
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Ia=Il+If 

Pl= la* Ia*Ra+Pk 

Pout  =V* II 

effi=Pout/(Pout+Pl) 
disp ( ef f i  ) 


Scilab  code  Exa  10.17  Probleml7 

V  =  250 
n  =  800 
Ra=0 . 15 
Rse  =0 . 1 
R=Ra+Rse 

Ial =120 
Eal=V-Ial*R 

Ia2  =60 
Ea2=V-Ia2*R 
n2=Ea2/Eal*n/0.7 
disp (n2) 


Scilab  code  Exa  10.18  Probleml8 


R  =  1 . 2 

V  =  220 

Ea35=V-35*R 
n35  =  (475  +  400) /2 

V  =  200 

E0=V-35* (R+2) 
n=n35*E0/Ea35 
disp (n) 
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Scilab  code  Exa  10.19  Problcml9 


V  =  250 
Ia=200 
Ra=0 . 22 

Ea=V- Ia*Ra 
Pm=Ea* la 
Prl =600 
Pmout=Pm-Prl 
n=1250 

wm  =  2*°/„pi*n/60 
Tl=Pmout / wm 
disp (T1 ) 

Rf =125 
Psh=V* V/Rf 
Pein=V* Ia+Psh 
effi=Pmout/Pein 
disp ( ef f i  ) 


Scilab  code  Exa  10.20  Problem20 

P=25000 
n=1600 
V  =  250 
If =1 . 5 
Rf =V/If 
disp (Rf ) 

Ra=0 . 1 
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V  =  220 
Ia=P/V 
Ea=V- Ia*Ra 
If =0 . 875 
disp(If) 

Rf =V/If 
disp (Rf ) 

Pdev=Ea* la 
disp ( Pdev ) 

Tdev  =  Pdev/2/°/0pi/ii*60 
disp ( Tdev ) 


Scilab  code  Exa  10.21  Problem21 


V  =  230 

Ea=210 

Ia=40 

Ra=(V-Ea) /la 
disp (Ra) 

Pdev=Ea* la 
n=1200 

wm  =  2*°/„pi*n/60 
Tdev=Pdev/ wm 
disp (Pdev ) 
disp (Tdev ) 

n=n*V/Ea 
disp (n) 


Scilab  code  Exa  10.22  Problcm22 
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V  =  400 
Ia  =  50 
n  =  500 
Ra=0 . 5 
Ea=V- Ia*Ra 

Tdev=Ea*Ia/2/%pi/n*60 

Iastart  =75 
Rs=V/Iastart-Ra 
disp (Rs ) 

Tstart=Tdev*( Iastart/ la) ~2 
disp (Tstart ) 

n2  =200 

Ea200=Ea*n2/ n 
Rs=(¥-Ea200)/Iastart-Ra 
disp (Rs ) 


Scilab  code  Exa  10.23  Problem23 

nl =1200 
Rfcrit=260/4 
n2=1600 
Vf actor=n2/ n2 
Vnl =395 
disp (Vnl ) 


Scilab  code  Exa  10.24  Problem24 


Ia=500 
Ra=0 . 05 
Vb  =  2 
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Va= Ia*Ra+Vb 
Vt  =330 
disp (Vt ) 


Scilab  code  Exa  10.25  Problem25 

V  =  240 
Ea  =  V 
n  =  800 
1  =  16 

Ke=V/n/I 
K  =  I/n 

n=  sqrt (V*2/Ke/K/ sqrt (2) ) 
di sp (n) 

Ia=  sqrt (2) *K*n 
disp ( la) 


Scilab  code  Exa  10.26  Problem26 
V  =  600 

Pout  =60000 
ef f i =0 . 85 
Pin=Pout / ef f i 
Il=Pin/V 
Rf =100 
If =V/Rf 
Ia= II - If 
Ra=0 . 16 
Ea=V- Ia*Ra 

n  =  900 
nO=n*V/Ea 
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14  sr  =  nO/n-l 

15  disp (nO) 

16  disp(sr) 

17 

18  TL  =  Pin-Pout 

19  Pcu  =  la*  Ia*Ra 

20  Psh  =  If  *  If  *Rf 

21  Prl  =  TL-Pcu-Psh 

22  disp(Prl) 
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Chapter  11 

Synchronous  Machine 


Scilab  code  Exa  11.1  Problcml 

P=1000000 
Vrated=6600 
If  =60 

Xs_un sat =4700/ sqrt (3) /98 
disp (Xs_unsat ) 

Xs_adj=Vrated/sqrt (3) /143 
disp (Xs_adj ) 

Ia=P/sqrt (3)/Vrated 
pf =0 . 8 

theta=-acos (pf) 

Ia=Ia*(exp(%i*theta) ) 

Ef=norm(Vrated+%i*Xs_adj  *  la*  sqrt (3) ) 
di sp ( Ef ) 

VR=Ef /Vrated-1 
disp (VR) 


Scilab  code  Exa  11.2  Problem2 
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P=45000 
V  =  440 
pf =0 . 8 

Ia  =  P/ sqrt (3) /V*exp(%i*acos  (pf)  ) 
Vt=V/ sqrt (3) 

Rs  =  0 . 2 
Xs=l . 8 

Ef  =  norm(Vt-Ia*(Rs  +  °/0i*Xs)  ) 

If =Ef /85 
Pein  =  pf  *P 
Rf  =35 

Fl=If *If *Rf 

Pin=Pein+Fl 

Pcu=norm(3*Ia*Ia*Rs) 

Pshaft=Pein-Pcu 

Prl=1500 

Pshaft_net=Pshaft-Prl 

effi=Pshaft_net/Pin 

disp(Pshaft_net) 

disp(If) 

disp ( ef f i  ) 


Scilab  code  Exa  11.3  Problcm3 

V=12500 
Xs  =8 


Vt=V/ sqrt (3) 
Ef  =Vt 

Ef =Vt*l .2 
Ia= ( Ef - Vt ) /Xs 
theta  =  70pi  / 2 
pf  =0 
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Pe  =0 

Qe=-sqrt (3) *V*Ia 
disp ( la) 
disp(Pe) 
disp(Qe) 
disp (pf ) 

Ef  =  Vt  *0 . 8 
Ia=(Vt-Ef ) /Xs 
theta  =  -°/„pi  / 2 
pf  =0 
Pe  =0 

Qe  =  sqrt  (3) *V*Ia 
disp ( la) 
di sp ( Pe ) 
disp(Qe) 
disp (pf ) 


Scilab  code  Exa  11.4  Problem4 

Pe =10000000 
V=12500 
Xs  =8 

Vt=V/ sqrt (3) 

Ef =Vt*l .2 

delta=asin(Pe/3/Vt/Ef*Xs) 
Ia=(Ef*exp(y„i*delta)-Vt)/yoi/Xs 
pf=real (Ia)/norm(Ia) 

Qe=-sqrt (3)*V*imag(Ia) 

disp (norm ( la) ) 
disp (pf ) 
di sp ( Qe ) 
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Scilab  code  Exa  11.5  Problem5 


P  e  i  n  =  1 0  0  0  0  0  0 
pf =0 . 9 
Xs  =3 . 24 

theta=acos (pf ) 

V=3300 

Ia  =  Pein/  sqrt  (3)  /  pf/V*exp  (°/0i*theta) 
Vt=V/ sqrt (3) 

Ef  =  norm(Vt-0/0i*Ia*Xs) 
Pemax=3*Vt*Ef/Xs 
Ia=(Vt  +  Ef*°/„i)/°/„i/Xs 
Qe=-sqrt (3) *V*imag(Ia) 

di sp ( Pemax ) 
disp (norm ( la) ) 
disp (real (Ia)/norm(Ia) ) 
di sp ( Qe ) 


Scilab  code  Exa  11.6  Problcm6 

Pe  =20000 
Pm  =  20000 
V  =  400 
pf =0 . 8 
Xs  =4 . 5 

Ia  =  Pm/sqrt  (3)  /V/pf*exp(°/0i*acos  (pf)  ) 

Vt=V/ sqrt (3) 

Ef  =Vt -%i  *  Ia*Xs  ///calculation  mistake  in  the 

book  at  this  step 

Efline=norm(Ef)*sqrt (3) 
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disp(Efline) 

disp(atan(iraag(Ef)/real(Ef))) 
delta  =  -°/„pi /2 

Efmin=Pe/3*Xs/Vt*sin(-delta) 
Ef minline=Ef min* sqrt (3) 
disp(Efminline) 
disp  (delta/ °/„pi  *180) 

Ia=(Vt  +  %i*Efmin)/°/„i/Xs 

disp (norm ( la) ) 

disp (real (Ia)/norm(Ia) ) 


Scilab  code  Exa  11.7  Problcm7 


f  =50 
ns  =  100 
P=110000 
pf =0 . 8 

p=120*f /ns 
disp (p) 
kVA=P/pf/1000 
di sp ( kVA ) 
kW=P/0 . 971/1000 
disp (kW) 

Tpm=kW *1000*60/2/ %pi/ns 
disp(Tpm) 


Scilab  code  Exa  11.8  Problcm8 

P  =  1 2  E6 
Q  =  6  E6 
V=22000 
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Xs  =8 

S  =  P  +  °/0i  *  Q 

theta=atan(C]/P) 

disp(theta/°/„pi*180) 

Ia=norm(S)/sqrt (3)/V 

Ef=V/sqrt(3)+°/0i*Xs*Ia*exp(-0/0i*theta) 

delta=atan(imag(Ef)/real(Ef)) 

disp  (delta/ °/„pi  *180) 

emf=norm(Ef)*sqrt (3) 
di sp ( emf ) 


Scilab  code  Exa  11.9  Problem9 

Isc=60 
If  =50 

Voc  =  15300 
Isc=60 

Xs_unsat=Voc/sqrt (3) / I  sc 
disp (Xs_unsat ) 

V=11000 
Isc=67 . 5 

Xs_sat=V/sqrt (3)/Isc 
disp (Xs_sat ) 

pf =0.85 

theta=acos (pf ) 

Ia  =  50*exp(-°/„i*theta) 
Vt=11000 
Vl=Vt / sqrt (3) 
Ef=Vl+%i*Xs_sat*Ia 
Efline=norm(Ef)*sqrt (3) 
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disp(Efline) 


Voc=Ef line 
VR=Voc/Vt-l 
disp (VR) 

If =57.5 
di sp ( If  ) 


Scilab  code  Exa  11.10  ProblcmlO 

P=10000 
V  =  400 
pf =0 . 8 
Xs  =  1 6 

theta=acos (pf ) 

Ia  =  P/sqrt  (3)/V*exp(-°/>i*theta) 

Vt=V/ sqrt (3) 

Ef  =  Vt  +  °/0i*Xs*Ia 
disp(norm(Ef)) 

disp(atan(imag(Ef)/real(Ef))*180/°/0pi) 

Ef  2  =  1 . 2*  norm ( Ef ) 

Pe=P*pf 

delta=asin(norm(Pe/3*Xs/Ef2/Vt)) 

Ef2  =  Ef2*exp(°/0i*delta) 

Ia=  (Ef  2 -Vt ) /°/0i/Xs  //calculation  mistake  in  the 

book  at  this  point 

disp (norm ( la) ) 
pf=real (Ia)/norm(Ia) 
disp (pf ) 

disp(acos (pf) *180/°/„pi) 
delta  =  70pi  / 2 

Pemax=norm(3*Ef*Vt/Xs*sin(delta)) 
disp ( Pemax ) 


182 


25 

26 

27 

28 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 


Ef =norm (Ef ) *%i 

Ia=(Ef-Vt)/°/0i/Xs 

disp (norm ( la) ) 

disp (real (Ia)/norm(Ia) ) 


Scilab  code  Exa  11.11  Problemll 


f  =50 

MechLoad= (8+0 . 5) *1000 

Pein=MechLoad 

Vt  =23 1 

Ef  =750/ sqrt (3) 

Xs  =  1 6 

delta=asin(Pein/3*Xs/Ef/Vt) 

Ef=Ef*exp(-%i*delta) 

Ia=(Vt-Ef  )/°/„i/Xs 

ns =120*50/4 

ws  =  ns/60*2*°/„pi 

Tdev=Pein/ ws 

Tshaft=8000/ws 

disp (norm ( la) ) 

disp (real (Ia)/norm(Ia) ) 

disp  (acos  (real  (la)  /  norm  (la)  )  *180/°/0pi) 

disp ( Tdev ) 

disp (Tshaf t ) 

Ef  =600/ sqrt  (3) 
delta=0 

Ia=(Vt-Ef)/%i/Xs 

kVAR  =  sqrt (3) *400* norm ( la) 

disp (norm ( la) ) 

disp (real (Ia)/norm(Ia) ) 

di sp ( kVAR ) 

C  =  norm(Ia)/Vt/2/°/0pi/f 
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Ef  =300/ sqrt (3) 

Ia=(Vt-Ef  )/°/„i/Xs 

kVAR  =  sqrt (3) *400* norm ( la) 

disp (norm ( la) ) 

disp (real (Ia)/norm(Ia) ) 

di sp ( kVAR ) 

L  =  Vt/norm(Ia)/2/°/0pi/f 
kVAR=6000 

Ia=kVAR/ sqrt (3) /400*%i 
Ef  =  Vt  -  °/0i  *  Xs  *  I  a 
disp(Ef*sqrt (3) ) 


Scilab  code  Exa  11.12  Probleml2 

P=1000000 
V=6600 
Xs  =25 
pf =0 . 8 

theta=  acos (pf ) 

Ia  =  P/sqrt  (3)/V*exp(-°/0i*theta) 

Vt=V/ sqrt (3) 

Ef  =  Vt  +  °/„i*Xs*Ia 

delta=asin(norm(Xs*Ia*cos(theta)/Ef)) 

Vt  =  norm(Ia)*Xs*sin(°/0pi-asin(pf)-delta)/sin(delta) 
Vl  =  ¥t  *  sqrt  (3) 
disp (VI ) 


Scilab  code  Exa  11.13  Problcml3 


Xs  =8 
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Ia=40 
V  =  600 

Vt=V/ sqrt (3) 
Ef  =V 


delta=acos ((l-(Xs*Ia/V) ~2+(Vt/V) "2) /2/Vt*V) 
theta=-atan ( ( cos (delta) -Vt/V) / sin (delta) ) 
pf  =  cos (theta ) 

Pout=sqrt (3) *V*pf*Ia 
di sp (Pout ) 

theta=- theta 
Pin=Pout 
di sp (Pout ) 


Scilab  code  Exa  11.14  Problcml4 


V  =  400 

Vt=V/ sqrt (3) 

Ia  =  50 
pf  =  l 
theta=0 
Xs  =  1 . 3 

Ef  =  Vt  -  %i  *  Xs  *  I  a 

delta=atan(imag(Ef)/real(Ef)) 

disp(delta*180/%pi) 

Pelec=sqrt (3)*V*pf*ia 

Pmech=Pelec 

pf =0 . 8 

theta=  acos  (0.8) 

Ia  =  Pelec/sqrt(3)/V/pf*exp(°/„i*theta) 
Ef2=Vt-%i*Xs*Ia 
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If=0.9*norm(Ef2)/norm(Ef) 
di sp ( If  ) 


Scilab  code  Exa  11.15  Probleml5 

Xspu=0 . 8 

P=1000000 

V=3000 

Xs=Xspu*V*V/P 

Vt  =  V/ sqrt  (3) 

Pm=750000 
Pe  =  Pm 
pf =0 .85 

theta=acos (pf ) 

Ia  =  Pm/0.8  /  sqrt  (3)/V*exp(-°/0i*theta) 


Ef  =  Vt  -  %i  *  Xs  *  I  a 
Efline=sqrt (3) *norm(Ef) 
disp(Efline) 

delta=%pi /2 
Pmax=3*Vt*norm(Ef)/Xs 
di sp ( Pmax ) 
ns =120*50/ 12 
ws  =  ns*2*°/0pi/60 

Tmax  =  Pmax/ws  ////// calculation  mistake  in  the 

book  at  this  point 

di  sp  ( Traax ) 

delta=90 

Efmin=Pm/3*Xs/Vt 
Efmin_line=sqrt (3) *Efmin 
disp(Efmin_line) 
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30  Ia=(Vt  +  %i*Efmin)/°/„i/Xs 

31  disp (norm ( la) ) 

32  pf =real ( la) /norm ( la) 

33  disp(pf) 
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Chapter  12 

Induction  Motor 


Scilab  code  Exa  12.1  Problcml 

V  =  400 
1=80 
pf =0.75 

Pin=sqrt (3) *V*I*pf 
Pshaf t=Pin*0 . 85 
disp (Pshaf t ) 
ns=1000 
s  =  0 . 04 
n=(l-s) *ns 
w  =  2*  °/0pi  *n/60 
Tshaft=Pshaft/w 
disp (Tshaf t ) 


Scilab  code  Exa  12.2  Problem2 


f  =50 

f 2=120/60 
s=f 2/f 
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ns=1500 
n=(l-s)*ns 
w  =  2*  °/0pi  *n/60 

T  =  1 00 

Pshaf  t  =T*w 
disp (Pshaf t ) 

Pm= ( T+7) *w 

Pcur=Pm*s/(l-s) 

disp ( Pcur ) 

Pin=Pm+Pcur+700 

disp(Pin) 

effi=Pshaft/Pin 

disp ( ef f i  ) 


Scilab  code  Exa  12.3  Problcm3 


ns =120*50/20 
sf 1=1 -292 . 5/ns 

R2  =0 .  12 
X2=l . 12 
smaxT=R2/X2 
di sp ( smaxT ) 

Tmax_by_Tf 1=0 .5/X2/(R2/sfl/((R2/sfl)“2+X2“2)) 
disp(Traax_by_Tfl) 


Scilab  code  Exa  12.5  Problcm5 


ns  =750 
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R2  =  0 . 5 
X2  =5 

smaxT=R2/X2 
n= ( 1 - smaxT ) *ns 
disp (n) 

k=R2/X2 

Tmax_by_Tstart=(l+k*k) /2/k 
disp(Tmax_by_Tstart) 

p  =  poly  (  [1  ,  -4 , 1]  ,”k”  coeff”) 
w=roots (p) 

R2=w*5 

R2ext  =  R2 -0 . 5 

R2total =18.66 
smaxT=R2total/X2 
R2ext  =  R2ext (2) 
disp(R2ext) 


Scilab  code  Exa  12.6  Problcm6 


Traax  =200 
Tstart  =80 

p=poly([l,-2*Tmax/Tstart,l]  ,  ”  smaxT”  ,  ”  coeff”  ) 
smaxT  =  roots (p) 
smaxT  =  sraaxT  (2) 
di sp ( smaxT ) 


p  =  poly  (  [1  ,  -4 , 1]  ,”w”  coeff”) 
w=roots (p) 
sf 1  =  smaxT /w  ( 1 ) 
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disp(sfl) 

ratio=sqrt (((smaxT/sfl) "  2  +  1 ) / (smaxT~2+l) ) 
disp(ratio) 

p  =  poly( [1 , -((smaxT/sfl) “  2  +  1 ) *  sf 1 , smaxT  ~  2]  , ” k”  , ” coeff 
”) 

k=roots (p) 
k=k (2) 
di sp (k) 


Scilab  code  Exa  12.7  Problem7 

V  =  440 
ns=1000 

ws=2*%pi*ns/60 
n  =  975 
s=l -n/ ns 

Z  =  1 . 06  +  0 .576/s  +  °/0i*(l.  68  +  0 .75) 

I2=V/ sqrt (3) /Z 

Im  =  -°/„i*V/sqrt  (3)  /44 . 2 

11= Im+ 12 

pf=real (Il)/norm(Il) 

Pin  =  sqrt (3)*V*norm(Il)*pf 

Pout  =  norm (3*12*12*0. 576*(l/s-l)) -415 

ef f i =Pout /Pin 

Tnet=Pout/ws/(l-s) 

di sp ( Tnet ) 
di sp ( Pin ) 
di sp ( norm  ( 1 1 )  ) 
disp (pf ) 
disp ( ef f i  ) 


191 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 


Scilab  code  Exa  12.8  Problem8 


V  =  400 

P0=1210 

10=8.7 

Ri = ( V/ sqrt (3) ) ~2/P0*3 
disp (Ri ) 

pf=P0/sqrt (3) /V/IO 
thetaO=acos (pf) 
Xm=Ri/tan(thetaO) 
disp (Xm) 

PI =6050 
VI =200 
11=47.8 
R=Pl/3/Il/Il 
RI =0 . 42 
disp (RI ) 

R2=R-R1 
disp (R2) 

Z=V1 / sqrt (3) /II 
X=sqrt (Z*Z-R*R) 
disp (X) 

ns  =750 
n  =  710 

ws  =  ns*2*°/0pi/60 
s=l -n/ ns 
Z  =  Rl+R2/s  +  °/0i*X 
I2=V/ sqrt (3) /Z 
I0=10*exp(-%i*theta0) 
11=10+12 
disp(norm(Il)) 
pf=real (Il)/norm(Il) 
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disp (pf ) 

T=norm(3*I2*I2*R2/s/ws) 
disp (T) 

Pout=T*ws*(l-s) 

Pin  =  sqrt (3)*V*norm(Il)*pf 
ef f i =Pout /Pin 
disp ( ef f i  ) 

s  =  l 

Z  =  R  +  %  i  *  X 
I2=V/ sqrt (3) /Z 
11=10+12 
disp (norm (II)) 

Tstart=norm(3*I2*I2*R2/s/ws) 
disp (Tstart ) 

smaxT=R2/ sqrt (R1 ~2+X~2) 

Z  =  R1+R2/ smaxT  +  °/„i  *X 
I2=V/sqrt (3)/norm(Z) 
Tmax=3*I2*I2*R2/ smaxT/ws 
disp ( Tmax ) 
di sp ( smaxT ) 


Scilab  code  Exa  12.9  Problem9 

Is_by_If 1=5 
sf 1=0 . 04 

Ts_by_Tfl=Is_by_Ifl ~ 2*  sf 1 
disp (” pu”  , T  s  _by _Tf 1 ) 


Scilab  code  Exa  12.10  ProblcmlO 
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R2  =  0 . 5 
X2  =5 

Tm_by_Tf 1=3 . 5 

Ts_by_Tfl=R2*Tm_by_Tfl*2*X2/(R2~2+X2~2) 

disp(Ts_by_Tfl) 

Ts_by_Tf 12=Ts_by_Tf 1/3 
disp(Ts_by_Tfl2) 

Ts_by_Tfl3=0.7*0.7*Ts_by_Tfl  //calculation 

mistake  in  the  book  at  this  point 
disp(Ts_by_Tfl3) 


Scilab  code  Exa  12.11  Problcmll 


si =1 -960/ 1000 
s2=l-800/1000 
R2ext=4/3*s2*0.25/sl-0.25 
disp(R2ext) 


Scilab  code  Exa  12.12  Probleml2 


smaxT=l-860/1000 
R2  =0 . 25 
X2=R2/ sraaxT 
Traax=180 
k=Tmax*X2 

s  =  0 . 045 
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T  =  2*  100 . 6/ (  (  R2/s )  ~ 2  +  X2  ~2)  *R2/ s  //some  mistake  in 

the  book  at  this  point.  (What  is  Vl/a) 

T  =  1 05 
disp (T) 

R2ext  =  X2 -R2 
di sp ( R2ext ) 


Scilab  code  Exa  12.13  Problcml3 
V  =  400 

Poutmech=5000 
Prl =285 

Pm=Poutmech+Prl 
s=l -1445/1500 
s2=l-900/1500 


p=poly([0,-(l-s)*V*V/Pm/s,l/s/s+l/s2/s2], ”R2” , ” coeff 
”) 

r=roots (p) 

R2  =  r  (1) 

X2=R2/s2 

ws  =2*  %pi *  1500/60 
Tmax=3/ws*V*V/3*0.5/X2 
disp ( Traax ) 


Scilab  code  Exa  12.14  Problcml4 

Pinout  =20000 
Pwf 1=1500 
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Pm=Pmout +Pwf 1 
disp (Pm) 
ns=1500 
n= 1440 
s=l -n/ns 
Pg=Pm/(l-s) 
disp (Pg) 
Pcu=s*Pg 
disp(Pcu) 


Scilab  code  Exa  12.15  Probleml5 

V 1 =400 
Rl =2 . 28*3/2 
Rlac  =  1 . 1*3 . 42 

Y0  =3 . 5/ sqrt (3) /VI 
Gi=(445/3) /Vl/Vl 
Bm=sqrt (Y0~2-Gi~2) 

Ri =1/ Gi 
Xm=l/Bm 

V2  =  200 

Z  =  V2  /  1 6 . 7*  sqrt  (3) 

R=2220/16 . 7/16.7 
X  =  sqrt (Z*Z-R*R) 

R2=R-Rlac 

n  =  935 
ns=1000 
s=l -n/ ns 

Z  =  Rlac  +  R2/s  +  °/„i*X 
I2=V1/Z 

Ii=Vl/Ri-'/.i*Vl/Xm 
I 1 = 12+ I i 
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Iline=norm(Il)*sqrt (3) 

pf=real (Il)/norm(Il) 

ws=2*%pi*ns/60 

T=norm(3*I2*I2*R2/s/ws) 

disp (T) 

disp(Iline) 

Pm=norm(3*I2*I2*R2*(l-s)/s) 

Pin=sqrt (3) *Vl*Iline*pf 

effi=Pm/Pin 

disp ( ef f i  ) 

V  =  400 
s  =  l 

Z  =  Rlac  +  R2/s  +  0/0i*X 

I2  =  V/sqrt  (3) /Z  ////calculation  mistake  in  the 

book  at  this  step 

Ts=norra (3*12*12* R2/ws) 
disp (Ts ) 


Scilab  code  Exa  12.16  Probleml6 
V  =  400 

Pinout  =5000 

Prl =285 

Pm  =  Pinout +  Prl 

s=l -1445/1500 
smaxT=l-900/1500 

R2=(l/s-l) *V*V/ Pm / (1/ s~2+l/smaxT~2) 

X2=R2/ smaxT 

ws=2*°/„pi*1500/60 

Tmax=V*V/ws/2/X2 

disp ( Tmax ) 
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Scilab  code  Exa  12.17  Probleml7 


V  =  440 
1=25 
pf =0.85 

Pin  =  sqrt (3)*V*I*pf 

ns=1500 

Pcore  =750 

Pcus  =950 

Pcur  =450 

Pwf 1=250 

Pg=Pin -Pcore -Pcus 
disp (Pg) 

Pm=Pg-Pcur 
disp (Pm) 

Pmout=Pm-Pwf 1 
di sp  ( Pinout ) 

effi=Pmout/Pin 
disp ( ef f i  ) 

s=Pcur /Pg 
n=(l-s) *ns 
disp (n) 


w  =  2*  °/0pi  *n/60 
Tdev=Pm/ w 
disp ( Tdev ) 
Tnet  =  Pmout /w 
disp (Tnet ) 
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Chapter  14 

Measurement  Techniques  and 
Electric  and  Electronic 
Instrumentation 


Scilab  code  Exa  14.1  Probleml 

Rvl=60/1000*10000 

Rv 2=120/ 1000*10000 

Rx  =  (Rv2-Rvl) *(1/(27.5/30*2-1)  -1) 

1=27 . 5/1000/600* (Rx  +  600)  /Rx 

Vact  =Rx* I 

disp (Vact ) 


Scilab  code  Exa  14.2  Problem2 


If sd=25/1000 
lx  =4 . 975 
Rx=Ifsd*5/Ix 
disp (Rx) 
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Rx=20/Ifsd-5 
disp (Rx) 


Scilab  code  Exa  14.3  Problcm3 


Rl=40/0 . 025-5 
R2  =  60/0 . 025-5-R1 
disp (R1 ) 
disp (R2) 


Scilab  code  Exa  14.4  Problem! 


W1 =500 
W2  =  -200 
P=W1+W2 

phi = at an ( ( W1 -W2) / ( W1+W2) ) 
pf=cos (phi) 
disp (pf ) 
disp (P) 


Scilab  code  Exa  14.5  Problem5 


E  =  4 

R1 =800 
R2  =80 
R3=1605 
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R4=160 
Rg  =  80 
s  =  8 


Eth=E* (R2/(R2+R4) -R1/(R1+R3) ) 
Rth=Rl*R3/ ( R1 +R3 ) +R2 * R4 / (R2+R4) 
Ig=Eth/ (Rth+Rg) 

D  =  s  * Ig  *  1 000000 
disp (D) 


Scilab  code  Exa  14.6  Problem6 

function  Zeq=parallel ( Z 1 , Z2) 
Zeq=Zl*Z2/ (Z1+Z2) 
endf unction 

f =1000 
w  =  2*  °/0pi  *f 

///AB 

R=2000 
C=0 . 045E-6 
Xc  =  l/°/„i/w/C 
Zl=parallel (R,Xc) 

//// BC 

R=1000 
C=0 . 45E-6 
Xc  =  l/°/„i/w/C 
Z2=R+Xc 


////AD 

C=0 . 4E-6 
Z3  =  l/°/0i/w/C 
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Z4=Z2*Z3/Z1 
di sp ( Z4 ) 


Scilab  code  Exa  14.7  Problcm7 


ksine=l . 11 
ksqr =1 

err=(ksine-ksqr)/ksqr*100 
disp(”  percent” , err) 


Scilab  code  Exa  14.8  Problem8 


f =50E6 
w  =  2*  °/0pi  *f 
theta  =  5*°/0pi/180 


C=1E-12 
R=l/ theta/ w/ C 
disp (R) 

C=10E-12 
R=l/ theta/ w/ C 
disp (R) 

C=100E-12 
R=l/ theta/ w / C 
disp (R) 
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Scilab  code  Exa  14.9  Problcm9 


f  =8E6 
w  =  2*  °/„pi  *f 
C=150E-12 
R  =  5 


Ql=l/w/C/R 
Rins  =0 . 1 

Q2=l/w/C/ (R+Rins) 

err=(Ql-Q2)/Ql*100 

disp(err) 


Scilab  code  Exa  14.10  ProblcmlO 

f =100E3 
w  =  2*  °/0pi  *f 
C=400E-12 

L=l/C/w/w 
disp (L) 


Scilab  code  Exa  14.11  Problcmll 


Resolution=l/10~4 
disp (Resolution) 
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disp (” 14.760” ) 


R_lV=l*Resolution 

di  sp  (”  0 . 5434  ”  )  ///since  any  digits  upto  4th  decimal 
can  be  shown 

R_10V  =  10*  Resolution 

disp(”  00.543”)  ///since  only  three  decimal  digits 
can  be  shown 


Scilab  code  Exa  14.12  Problcml2 


Vp  =  10 
Vref  =5 
R=100000 
C=0 . 22E-6 
T=Vp*R*C/Vref 
disp (T) 


Scilab  code  Exa  14.13  Problcml3 


V=3 .217 
bits=4 
base  =5 
a  =  0 


for  i=l : bits 

a=floor(V/base*(2~i)) 
disp (a) ; 

V=V-a*base/ ( 2  ~ i ) 
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end 


Scilab  code  Exa  14.14  Probleml4 


reading=5 

x=0 . 005/ 100* reading  / / 0.005%  of  reading 

err  =  x  +  lE-6  //final  digit 

disp (” seconds”  ,  err) 

reading=500 

x=0 . 005/ 100* reading  / / 0.005%  of  reading 

err  =  x  +  l  //final  digit 

disp(”  seconds”  ,  err) 


Scilab  code  Exa  14.15  Problcml5 


V  =  5 . 9 
Rx  =  0 
Rm=2000 
R1 =49000 
If sd=100E-6 
Vara=Rm*Ifsd 
1= (V- Vara) /R1 
Ish=I - If sd 
R2=Vam/Ish 
disp (R2) 


1=0 . 6*1 


Req=V/I 

Rx=Req-Rl -Rm*R2/ (R2+Rm) 
disp (Rx) 
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Scilab  code  Exa  14.17  Probleml7 


f  =50 

w  =  2*  °/„pi  *f 

Z2=1000 
Z3  =  1 6800 

C=0 . 38E-6 
Xc  =  l/°/„i/w/C 
Z4  =  833  +  Xc 

Z 1  =  Z2  *  Z3  /  Z4 
Rx=real ( Z 1 ) 
Lx= imag ( Z1 ) /w 

disp (Rx) 
di sp ( Lx ) 


Scilab  code  Exa  14.18  Probleml8 


f  =50 

w  =  2*  °/0pi  *f 
T  =  1  /  f 

Iav  =  l/ (T/2) ^integrate (  ’  sin (w* t)’,’t’,0,T/2) 
disp ( lav) 
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Chapter  15 
Power  Systems 


Scilab  code  Exa  15.1  Problcml 

VI =250 
V2  =480 

Vol2_by_Voll=Vl/V2 

sav=(l-Vol2_by_Voll ) *100 
di sp ( sav ) 


Scilab  code  Exa  15.2  Problem2 

P  =  5E6 
pf =0 .85 
V=33000 
1=50000 
rho  =3E-8 
Pt=P*pf 
Pl=Pt*0 . 1 
I  =  P/V 

Al=2*I*I*rho*l/Pl 
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Voll=2*l*Al 
disp (Voll ) 

Il=P/sqrt (3) /V 
A2=3*Il*Il*rho*l/Pl 
Vol2=3*l*A2 
disp ( Vol2 ) 


Scilab  code  Exa  15.3  Problcm3 


f  =50 

w  =  2*  °/0pi  *f 
1=0 . 8 
V  =  220 
P  =  75 

phi=acos (P/V/I) 
phi_new  =  acos  (0.9) 

Ic=I*cos (phi) *(tan(phi) -tan(phi_new) ) 

C=Ic/V/w 

disp (C) 

phi_new=acos (1) 

Ic=I*cos (phi) * ( tan ( phi ) -tan(phi_new) ) 

C=Ic/V/w 

disp (C) 


Scilab  code  Exa  15.4  Problem! 

Cond_cost  =  100 
charge  =60 

phi2=asin (0 . l*Cond_cost/ charge) 
pf=cos (phi2) 
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disp (pf ) 


Scilab  code  Exa  15.5  Problem5 


0c =400000 
pf 1=0.8 

phi 1  =  acos ( pf 1 ) 
ab=0c/cos (phil)*sin(phil) 
pf 2=0 . 25 
phi3=acos (pf 2 ) 
pf 2=0 . 484 

gammaa= (ab-pf2*0c) / (pf2*cos (phi3)+sin(phi3) ) 
di sp ( gammaa ) 


Scilab  code  Exa  15.6  Problcm6 


f  =50 

w  =  2*  °/0pi  *f 
P  =  2E6 
V=11000 
pf =0 . 8 

phi=acos (pf ) 

Xl  =  10 

IR=P/sqrt (3)/V/pf 
Vr=V/ sqrt (3) 
Vs=Vr+IR*Xl*sin(phi) 
Vsll=Vs*sqrt (3) 
disp ( Vsll ) 
VR=Vsll/V-l 
disp (VR) 
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Pf  =  1 

disp (pf ) 
Qc=P*tan (phi ) 
C=Qc/V/V/w 
disp (C) 


Scilab  code  Exa  15.7  Problem7 


f  =50 

w  =  2*  °/0pi  *f 
V=33000 
Vr=V/ sqrt (3) 

P  =  24  E6 / 3 
pf =0 . 8 

phi=acos (pf ) 

Ia=P/Vr / pf 
R1  =4 
XI  =20 

Vs=Vr+Ia*(Xl*sin(phi)+Rl*cos (phi) ) 

Vsll=sqrt (3) *Vs 
VR=Vsll/V-l 
disp ( Vsll ) 

Ia= la*  exp ( -%i *phi ) 
disp (norm ( la) ) 

phil=atan(-Rl/Xl) 
pf=cos (phil) 

Ial=P/Vr/ pf 

Ial  =  Ial  *  exp  ( -°/„i  *phi  1 )  ////calculation  mistake 

in  the  book  at  this  step 

Ic=Ial -la 
C=norm(Ic/w/Vr) 
disp (C) 
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27  LL1 =norm ( la* Ia*Rl ) 

28  effil=P/ (P+LL1) 

29  LL2  =  norm ( Ial * Ial *R1 ) 

30  effi2=P/ (P+LL2) 

31  disp(effil) 

32  disp(effi2) 
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